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ABSTRACT 


The ability to produce an accurate estimate of the expected amount of damage that 
will occur to a targeted region from the fragments generated by intercepting a Theater/ 
Tactical Ballistic Missile (TBM) in flight is necessary for determining whether the TBM 
should or should not be intercepted. For a sea-based Theater Missile Defense (TMD) 
system that is capable of predicting the remaining portion of the TBM's trajectory along 
with its impact point from its current and past flight parameters, possessing this ability 
may prevent unnecessary intercepts, and thus unnecessary damage. This thesis proposes a 
simplistic simulation model that produces expected damage estimates for various intercept 
ranges. Damage is evaluated in three areas: number of fragments to impact the target 
region, total mass to impact the target region, and total kinetic energy to impact the target 
region. The values generated by this model can aid in the decision of whether to intercept 
the TBM or not by comparing them against the expected damage values caused by the 
unintercepted TBM which are calculated by the sea-based TMD system. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may 
not have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic errors, 
they cannot be considered validated. Any application of these programs without 
verification is at the risk of the user. 
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EXECUTIVE SUMMARY 


PROBLEM: A model for estimating the expected collateral damage to a target region 
resulting from the intercept of a Theater/Tactical Ballistic Missile (TBM) is necessary for 
the decision maker to determine if intercepting the TBM is the proper course of action 
and is the focus of this thesis. 

With the increased threat of theater/tactical ballistic missile attacks against 
defenseless regions, and the propensity of some nations to use ballistic missiles as a 
bullying tactic, theater missile defense (TMD) has become a high-priority issue for the 
Department of Defense. Presently, TBM's pose a great threat to U.S. military operations 
because almost every conceivable future conflict involving U.S. military deployments 
would involve the risk of TBM attacks against U.S. forces. The only existing deployable 
TMD system which intercepts a TBM in flight. Patriot, an Army air defense system 
adapted for the TMD role, has suffered severe criticism for its performance in Desert 
Storm. At first, Patriot was praised as an extremely effective interceptor of SCUD 
missiles launched from Iraq at Israel and Saudi Arabia. Unfortunately, post-war analysis 
showed that the Patriot missile was much less effective than was initially claimed. Patriot 
is now being modified to enhance its capabilities against TBM's. Another land-based TMD 
system under development is the new Theater High Altitude Area Defense system 
(THAAD). The integration of THAAD and an improved Patriot promise to provide 
outstanding defense against future ballistic missiles. 

In addition to these systems, the Ballistic Missile Defense Organization (BMDO) 
has directed the Navy to adopt this two-tiered approach of exo-atmospheric and 
endo-atmospheric defense and develop its own TMD systems to be integrated with the 
Army TMD systems in the future. The weapon chosen for the lower-tier, or endo- 
atmospheric, role is SM-2 Block IV A surface-to-air missile. One concern with the future 
testing of this missile is the unknown magnitude of the collateral damage that results from 
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a TBM intercept. Therefore, a model is proposed for estimating the expected damage to 
a target region and its surrounding areas. 

SOLUTION METHODOLOGY : This thesis addresses the question of when to 
intercept a TBM and when not to intercept it. The model proposed is a simplistic 
analytical attempt at estimating the distribution of the expected damage that occurs to a 
target region and its surrounding areas from intercepting a TBM at predetermined points 
on its flight path. The estimate can then be compared against the amount of damage that 
would have occurred if the TBM were allowed to continue without intercept, assuming 
the TMD system in question has this capability. 

Utilizing MODSIM, a simulation, which replicated the behavior of the fragments 
generated after an intercept, was developed. Because of the lack of available data 
concerning the true fragmentation of a TBM upon intercept, many simplifying assumptions 
were made. In order to account for the uncertainty in the TBM's impact point, if not 
intercepted, two parameters of the TBM's initial conditions were allowed to vary in each 
run: vertical launch angle and off-target launch angle. Replications of the simulation were 
conducted for various intercept points and pairs of launch angles so as to generate data on 
the number of hits, total mass and total kinetic energy impacting the target region. The 
distributions of the number of fragment hits, the total mass impacting the target region and 
the total kinetic energy impacting the target region are examined to determine expected 
values. 

SOLUTION: The modeling and simulation results presented are speculative and 
exploratory. They illustrate the effects of various physical influences upon missile 
collateral damage to a target in a preliminary way. The conclusion drawn from this thesis 
is that knowing the expected damage to occur to a target region if a TBM is intercepted is 
helpful in making a decision whether to intercept it or not. This model is a simple 
analytical attempt to produce a tentative estimate. At this time, very little is known 
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concerning the dispersion of fragments generated by an intercept. Because of this lack of 
information, a crude estimate is the only estimate available. When the assumptions made 
in this thesis are determined with certainty, a more detailed version of the model, 
replicating the exact flight paths and impact points, can be developed. 
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I. INTRODUCTION 


With the increased threat of theater/tactical ballistic missile (TBM) attacks against 
defenseless regions, and the propensity of some nations to use ballistic missiles as a 
bullying tactic, theater missile defense (TMD) has become a high-priority issue for the 
Department of Defense. Presently, TBM's pose a great threat to U S. military operations 
because no truly effective ballistic missile defense systems exist. Moreover, almost every 
possible future conflict involving U.S. military deployments would also involve the risk of 
TBM attacks against U.S. forces. The most familiar TBM, encountered extensively 
during the Gulf War, is the SCUD, Variant B. Many countries have these missiles in their 
arsenals including such unstable nations as: North Korea, Iraq, Vietnam, Yemen, Libya, 
Afghanistan, and Syria [Ref. 1], 

The only existing TMD system which intercepts a TBM in flight, the Patriot, an 
Army air defense system adapted for the TMD role, has suffered severe criticism for its 
performance in Desert Storm. At first, Patriot was praised as an extremely effective 
interceptor of SCUD missiles launched from Iraq at Israel and Saudi Arabia. 

Unfortunately, post-war analysis showed that the Patriot missile was much less effective 
than was initially claimed. In fact, the majority of SCUDs thought to be intercepted had 
spontaneously broken into segments upon re-entry as a result of faulty manufacturing. 
Patriot is now being modified to enhance its capabilities against TBM's. Another 
land-based TMD system under development is the new Theater High Altitude Area 
Defense system (THAAD). The integration of THAAD and an improved Patriot promise 
to provide outstanding defense against future ballistic missiles. However, THAAD 
development is being decelerated because of budget cuts. The drawback of both Patriot 
and THAAD, when the latter is completed, is the enormous lead time required to deploy 
and set up these systems. Future TBM attacks on areas of U.S. interest may not be 
known weeks in advance; time to deploy these defense systems may not be available. Also, 
permission must be granted from foreign leaders to deploy land-based TMD systems. The 
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leaders may be hesitant even though these systems are there to save their territory because 
deployment of such systems on their land may escalate hostilities. In addition to the 
deployment lead time is the enormous expense of airlifting the systems to the threat areas. 
In summary, such drawbacks provide strong motivation for the development of a 
sea-based ballistic missile defense system. 

As mentioned previously, plans and initial steps necessary for TBM attacks may 
not be known long enough (e g., weeks) in advance of attack to allow for deployment of 
Patriot batteries. A timely response option is the use of a mobile sea-based ballistic missile 
defense system. The Navy is in the process of developing a version of this system by 
upgrading the Aegis Weapon System and the SM-2 Block IV standard missile in order to 
effectively detect, track, engage, and intercept TBM's. The new standard missile is to be 
an improvement of the existing vertical-launch Terrier/standard surface-to-air missile. The 
latter requires modifications of the warhead, fuze, and seeker to improve its capability to 
defend against TBM's. The new version would build upon already-developed engineering 
and design efforts used on Patriot, but there are drawbacks to using this technology [Ref. 
1], One of these is the explosive warhead. This warhead may cause total destruction with 
a near-miss intercept, or it may not. If the near-miss intercept does not totally destroy the 
TBM, debris from the TBM could impact the intended target or surrounding areas. 

Collateral damage is the term for non-target damage caused by the debris. For 
example, if a TBM is targeted to impact the center of a coastal city and U.S. forces are 
conducting an amphibious assault, intercept of the TBM might stop it from directly 
impacting its target, but the fragments generated by the intercept might fall on the U.S. 
forces and interfere with the assault operation. When contemplating an intercept, the 
pattern of the debris generated from the intercept must be estimated in order to determine 
the amount of damage that will occur if the ballistic missile is intercepted and fragments 
continue on course. According to Israeli reports after Desert Storm, the 11 SCUDS 
intercepted by Patriot batteries caused more damage than did the 13 SCUDS that directly 
hit targets within Israel. This raises an interesting question: with prior knowledge as to 
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the amount of damage, both primary and collateral, that will occur from intercepting a 
TBM. is intercepting the TBM always the best strategy? If the TBM is allowed to 
continue on course without intercept, it is possible that the amount of damage resulting 
from the single impact point may be less than that of a shower of fragments impacting the 
target area. Of course, there is also a chance that the TBM will miss the target area 
completely. Because the dispersion of fragments resulting from an intercept and the size 
of each fragment are unknown, an accurate estimate of the debris pattern is almost 
impossible. This dispersion uncertainty may be caused by numerous aspects including 
angle of impact by interceptor, angle of descent of ballistic missile, impact point on 
ballistic missile, distance between interceptor and ballistic missile at intercept. 

Purpose of this Thesis: Where, and If, to Intercept? 

In order to address the question as to when to intercept a ballistic missile and when 
not to intercept it, a simulation model is developed to estimate the distribution of the 
damage that occurs from intercepting a TBM at predetermined points on its trajectory. 

This damage is then compared to the amount of damage that would have occurred if the 
TBM were allowed to continue without intercept. Development of the model requires a 
number of simplifying assumptions to be made. The first, and possibly most credible, set of 
assumptions deals with the characteristics of the TBM, including launch acceleration, 
motor burn time, missile length, missile diameter and missile mass. The second, more 
speculative, set of assumptions concerns the generation and motion of fragments, 
including number of fragments, mass of each fragment, cross-sectional area of each 
fragment, and the velocity of each fragment after intercept. These characteristics are 
generated from various probability distributions judged to be appropriate. The final 
assumption is target region designation. For this study the target region is a coastal region 
that is 2 km X 2 km. This region is then broken up into 16 blocks that are 0.5 km X 0.5 
km each. Actual targets that are subregions of this region can be studied for damage. 
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The preceding assumptions remain the same for each run of the model. However, 
two parameters are allowed to vary in each run: vertical launch angle and off-target 
launch angle. For each pair of angles, while holding all other parameters constant, a 
different TBM flight path is generated. From these flight paths an impact point could be 
calculated. Although, having precise knowledge of these two angles is not realistic, 
determining the impact point of a TBM from its flight path, conditional on the angles, is 
realistic For simulation purposes the two angles are taken to be given and known in order 
to calculate the impact point. Realistic expected impact damage can be estimated by 
averaging over any angle distribution of interest. 

The model is intended to represent a worst-case situation: the intercept results in a 
near-miss warhead explosion that creates missile fragments which continue traveling 
towards the target area with minimal loss of mass and minimal change in velocity. The 
TBM's flight path is modeled by using suitable equations of motion in three dimensions. 

At pre-determined points on the flight path, the hypothetical intercept points, the missile 
breaks up into a random number of fragments resulting from the near-miss intercept. Each 
fragment has the velocity of the TBM at intercept plus a random "kick" in three 
dimensions to simulate breaking apart. Each of the fragments continues towards the 
target region on a path determined by its mass and cross-sectional area. This model treats 
the cross-sectional areas of the fragments as functions of their respective masses, i.e., the 
larger the mass the larger the cross-sectional area. 

The objective of this modeling project is to estimate the expected damage that 
occurs when the TBM is intercepted at various distances from the target area, given the 
vertical and off-angle launch angles. Certain intercept points generate fragments which 
impact the target area. The amount of damage that occurs from a fragment impact is 
estimated using a damage function that takes into account fragment mass and velocity. 

The damage generated from different intercept points is then compared to the damage 
generated by letting the missile through the defenses with hope of a miss. Replications of 
this simulation are conducted for various intercept points and pairs of launch angles so as 
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to generate data on the numbers of hits, total mass and total kinetic energy impacting the 
target area. The distributions of the number of fragment hits, the total mass impacting the 
target area and the total kinetic energy impacting the target area will be examined to 
determine expected values and summaries of the effect of intercepting missiles at various 
locations along their trajectories. The ballistic missile modeled in these simulations are 
based on the characteristics of the SCUD variants used in Desert Storm, i.e., they will 
have a fairly large circular error probable (CEP) which is represented by the varying 
angles 

The modeling and simulation results presented are speculative and exploratory. 
They illustrate the effects of various physical influences upon missile collateral damage to 
a target in a preliminary way. The model equations and software, which are provided 
herein, can be modified or extended to test sensitivity to specific parameter values. With 
slightly more effort the actual physical relationships that have been used can, if considered 
necessary', be modified in more reasonable and defensible ways. 
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II. NATURE OF THE PROBLEM 


A. PROBLEM DESCRIPTION 

As demonstrated in the Gulf War, Iraq used it's A1 Hussein TBMs, i.e. modified 
SCUD B's, against civilian targets, most of which were populated cities. The circular error 
probable (CEP) of this particular TBM is fairly large, being about 1000 m, allowing a 
chance of missing the target altogether and landing in some insignificant area, e.g., the 
sea. Actually, this happened quite often during the Gulf War. Of the 88 TBMs launched 
by Iraq towards Israel and Saudi Arabia, 35 missed their targets completely, landing either 
in the sea or in the desert. Future sea-based TMD systems will be able to estimate, with a 
level of certainty, the impact point of a TBM by calculating the remaining portion of its 
trajectory according to the TBM's current and past parameters. However, these systems 
do not rely solely on the estimate in determining what action to take, since the estimate is 
not exact If the impact point is judged likely to be close to the target area, one of two 
courses of action must be taken: intercept the missile and risk collateral damage by the 
missile fragments, or do not intercept the missile and hope nothing significant is located 
near the impact point. Therefore, the problem to be addressed in this section is to 
determine whether a sea-based TMD system, assigned to protect a targeted region, e.g., a 
coastal city, should intercept a TBM in flight, creating missile fragments that may still 
impact the target area, or should refrain from interception , allowing the TBM to continue 
towards the intended target area, having assessed that it will miss. 

To address this problem, a computer simulation utilizing MODSIM, a simulation 
software system, is designed to determine the amount of damage that occurs from 
intercepting a TBM in flight. In this simulation a TBM modeled after the A1 Hussein, but 
with generic numbers substituted for some of its characteristics to preserve the 
unclassified status of this project, is launched toward a predefined target area. Parameters 
of the TBM are first chosen to allow an impact point at the exact center of the target area. 
Points along the TBM's flight path are selected as intercept points, at which the missile 
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separates into a random number of fragments that continue towards the target. The final 
position, or impact point, of each of the fragments generated is examined to determine if 
the fragment landed inside the targeted area, and where therein. Summing the results for 
each fragment totals the amount of damage to occur for a given intercept point. 
Conducting numerous replications of the model for each intercept point yields the average 
damage to occur to the target area. These results are further analyzed using statistical 
techniques to determine the expected damage to occur for various intercept points. 

Since this model is deterministic, the unintercepted TBM, under the assumption of 
a flat, non-rotating earth, would have the same impact point, i.e., the center of the target 
area, for each simulation run unless one or more parameters are allowed to vary. By 
changing parameters on each simulation run, the TBM's impact point, if allowed to 
continue towards the target without intercept, also changes. The parameters chosen to 
vary in this simulation are the vertical and off-target launch angles as depicted in Figure 1. 
By introducing a small degree of error in either angle, the TBM impacts vary from the aim 
point. As mentioned previously, a sea-based TMD system will be able to estimate, fairly 
accurately, this final hypothetical impact point of the TBM by constructing the remaining 
trajectory using its current and past parameters: course, speed, angle of descent, altitude, 
and others. Hence, a decision can be made either to intercept the TBM, and where along 
its path, or to allow it to continue on course by comparing the damage caused by a TBM 
impact at the point estimated by the sea-based TMD system against the calculated damage 
caused by fragments impacting the target area. 
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Figure 1. 3-Dimensional Plot of the TBM Launch Angles 

t = off-target launch angle 
v = vertical launch angle 


B. PROBLEM ASSUMPTIONS 

To preserve the unclassified status of this thesis and to allow the expressions 
utilized in developing the TBM trajectories to be closed form, the following assumptions 
are made: 

- The sea-based TMD system is stationed somewhere between the launch point 
and target where all potential intercept points utilizing an SM-2 Block IV A Standard 
missile are feasible. The position of the ship will not coincide with the debris pattern. 
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- Initial detection of the TBM is via space-based assets which provide a cue to the 
sea-based TMD weapons system. The sea-based TMD system picks up the TBM track as 
soon as it enters the detection envelope, giving it enough time to calculate an impact 
point. Error in that estimate will realistically occur, but is not considered in what follows. 

- The attack on the target consists of only one TBM. 

- The TBM experiences constant acceleration and straight-line motion during 
boost phase with air resistance having no effect. 

- The TBM follows a normal ballistic path, executing no maneuvers to avoid an 
intercepting missile. 

- The TBM remains intact throughout the flight and the mass remains essentially 
constant during boost phase even though fuel is being expended. At the completion of 
boost phase the TBM's mass is decreased by one seventh of its original mass to 
compensate for the spent fuel. 

- The intercept simulates a near-miss explosion of the SM-2's warhead which 
separates the TBM into several fragments, the number and sizes are simulated from a 
combination of plausible probability distributions; these latter can be modified as desired 
by an analyst. 

- There is minimal loss of mass and minimal change in velocity resulting from the 
intercept. 

C. MODEL OBJECTIVE 

If it is estimated that the TBM will impact in the vicinity of the target, intercepting 
it towards the end of its flight path at low altitude may generate fragments that still impact 
the target region and cause collateral damage. Intercepting it earlier in its flight path 
obviously presents a smaller chance of collateral damage. With this knowledge, the ideal 
intercept is the farthest distance possible from the target where intercept is feasible. 
However, at this distance the impact point of the TBM estimated by the sea-based TMD 
system is fairly inaccurate, and an unnecessary intercept may occur, i.e., the TBM would 
have fallen short of its target and landed in the sea, requiring no intercept. Before 
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planning an engagement, knowing the damage that results from different intercept points 
allows the sea-based TMD system to decide whether an intercept is necessary. Therefore, 
the objective of this model is, for particular intercept points, to estimate the expected 
values for the number of fragment hits , the amount of mass to impact the subregions of 
the general target region, and the total kinetic energy to impact subregions of the general 
target region in order to decide whether the damage sustained from the fragments will be 
greater than that of the TBM itself. The above measures of TBM effect are actually 
tabulated for subregions of the general target region. Note that the simulation allows 
evaluation of measures of target (sub)region effects other than expected values: it is, for 
example, possible to estimate the probability that the number of fragments that impact in a 
particular subregion exceeds a prescribed level, or the probability that the total kinetic 
energy delivered to a particular subregion exceeds a given level. This capability of the 
model is appropriate when hardened military targets are in the area, and their vulnerability 
to damage is low unless a specific threshold is exceeded. 
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III. MODEL 


The purpose of this thesis is to develop a simulation model that produces a 
tentative estimate of the expected damage that occurs within a designated target region 
caused by the fragments generated from intercepting a TBM. To accomplish this task the 
model may be considered as having two parts. The first part represents launching a TBM 
towards a designated target region located some known distance from a selected launch 
point. The TBM is tracked from launch through its trajectory until it reaches the intercept 
point chosen by the user. The effect of tracking the TBM is simulated by formulating and 
solving three-dimensional equations of motion, utilizing the TBM's mass and 
cross-sectional area and physical laws. 

The second part of the model begins with the TBM intercept, which is treated as a 
near-miss intercept that only breaks the TBM into pieces rather than destroying it. The 
reason behind this is to have the model represent a worst-case situation. If minimal mass 
and velocity are lost as a result of an intercept there stands a greater possibility of the 
fragments reaching the target region. At intercept the TBM breaks up into a random 
number of fragments. Each fragment generated continues towards the target region on its 
own flight path. These fragment flight paths are calculated using the same 
three-dimensional equations of motion as describe the TBM's flight path, but this time 
depending upon each fragment's mass and cross-sectional area. Of concern in this model 
is the extent of damage to the target region caused by these fragments. Therefore, the 
model determines which fragments impact the designated target region by comparing the 
final positions of each fragment to the target location. If it is determined that a fragment 
has impacted the target region, its mass and kinetic energy are recorded and subsequently 
totaled with the masses and kinetic energies of other impacting fragments. However, if it 
is determined that a fragment does not impact the target region, it is considered a miss, 
and the mass and kinetic energy are recorded for that region surrounding the target where 
the miss occurred. 
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Each run of the present model is analyzed by statistical techniques to estimate the 
expected number of fragments to impact various subregions of the target, and the 
expected total mass and expected total kinetic energy to impact these subregions. 
Conducting 2000 replications of the model for each intercept point enabled Monte Carlo 
analysis to be implemented to produce a fairly accurate estimate of these expected values. 
This many replications were required to keep the standard error of the estimate at a 
desired low level. The following section describes the model in depth. 

A. SIMULATION DESCRIPTION 

1. Initial Conditions 

As mentioned previously, the TBM utilized for this thesis is being modeled after 
Iraq's version of a modified SCUD, the A1 Hussein. For classification purposes, only four 
actual characteristics of the A1 Hussein are used in this simulation. The remaining 
characteristics required for the simulation have been calculated based on assumptions 
made concerning the distance to the target region. Another assumption concerning the 
flight of the TBM was that the calculations were based on a flat, non-rotating earth. For 
this thesis, the only portion of the TBM's trajectory of interest is immediately pre-intercept 
to impact. Therefore, the manner in which the TBM's flight path is calculated is artificially 
simplistic and thus is performed as a means of obtaining the trajectory portion of interest. 
The characteristics of the TBM modeled in this thesis, along with the calculations 
necessary to compute the simplistic TBM flight path are contained in Appendix A. 

2. Intercept and Fragment Generation 

In order to intercept a TBM, an actual interceptor of some kind is required. 
However, since the objective of this model is to determine the extent of damage to a target 
region caused by TBM fragments, given that an intercept has occurred, modeling of the 
interceptor is omitted. Therefore, to simulate the effect of intercept, the TBM breaks 
(spontaneously) apart into a random number of fragments, each experiencing a minimal 
change in velocity from that of the TBM because most of the energy created by the 
intercept is consumed fragmenting the TBM. Also, the accumulated mass of all the 


14 



fragments generated is taken to be close to that of the TBM. Lacking any data to validate, 
these assumptions are as good as possible. Accepting these two conditions, minimal 
changes in velocity and mass, it is plausible that the intercept represents a worst-case 
situation. 

Determining the exact number of fragments, and the size of each fragment, that 
result from intercepting a TBM is impossible because of numerous factors including angle 
of impact by interceptor, angle of descent of ballistic missile, impact point on ballistic 
missile, distance between interceptor and ballistic missile at intercept, and strength of 
TBM materials. Therefore, estimates of the random number of fragments and of the 
individual fragment characteristics, i.e., mass, cross-sectional area, and initial velocity, are 
introduced to complete the model. 

The actual realized numbers of fragments and the individual fragment masses that 
occur in one simulation replication are developed using a two-step operation. In the first 
step, a random number is drawn from the Poisson probability distribution. The mean 
number of fragments desired for this model, which is requested from the user at the start 
of the simulation, is used in generating this Poisson random number. This number is then 
divided into the TBM's mass to get the average mass of a fragment. The second step 
involves a loop which generates fragments one at a time. This loop is dependent solely on 
the mass of the fragments and continues until the accumulated mass of the fragments is 
greater than that of the TBM. Each time through the loop a fragment mass is generated 
by multiplying the average fragment mass produced in the first step by an Exponential 
random variable with a mean of one. This process simulates the variability between 
fragment masses. Each fragment mass is then added to the accumulated total. So long as 
the total is less than the TBM's mass, a fragment is generated. When the accumulated 
total becomes greater than the TBM's mass, fragment generation ceases and the number of 
fragments is determined. Of course modifications of the above procedure can be made in 
the simulation code: both Poisson and Exponential distributions are provisional illustrative 
choices. 
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The remaining characteristics, cross-sectional area and initial velocity, of each 
fragment are also calculated during each run through the loop. Three assumptions were 
made in order to specify the cross-sectional area of a fragment. First, since the exact 
shape of each fragment generated is unknown, this model proposes that "the" 
cross-sectional area is a function of the fragment's mass. Second, it assumes that each 
fragment, at the instant of intercept, takes on the (approximate) shape of a sphere. The 
last assumption is that the mass density of the fragment equals that of the TBM. With 
these assumptions, each fragment's cross-sectional area is calculated using simple 
geometric equations: 


Volume = ^ ass 

Density 

(18) 

Area = n ■ radius 2 

(19) 


Once again, because of the lack of knowledge concerning the way in which fragments are 
ge ■ mated after an intercept of a TBM, a presumption as to how the model should handle 
fragment velocities was necessary. The initial velocity of each fragment is fairly easy to 
calculate, assuming the fragments are dispersed in a Gaussian manner, i.e., the TBM 
breaks up in a conical fashion, and that there is negligible loss of velocity. To accomplish 
f his Gaussian dispersion, a random "kick" was applied to the three components of the 
aBM's velocity at intercept. The amount of this "kick" was determined by a Normal 
random variable with a mean of zero and standard deviation of ten meters/second. This 
allows each fragment to have a different initial velocity and, when combined with different 
masses and cross-sectional areas, a different flight path. 

3. Fragment Flight Paths 

The fragment flight paths are calculated in the same manner as is that of the TBM. 
A loop is initiated at intercept for each fragment with its initial conditions estimated under 
the above assumption. At each pass through the loop the fragment's velocity is calculated 
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using equations (12) - (14) from Appendix A. The fragments experience a significant 
reduction in their velocities. This change is the result of greater air resistance being 
applied to the fragments mainly because of the reduction in their respective masses. 
Fragments with larger masses experience smaller reductions in velocity than do those with 
smaller masses. The positions of the fragments along their paths are then updated using 
equations (15) - (17) from Appendix A. Prior to completing each pass through the loop 
an extra step, not used in determining the flight path of the TBM, is implemented. The 
extra step concerns the cross-sectional area. In calculating the TBM's flight path, which in 
this model is stable, the cross-sectional area remained constant. However, it is believed 
that the flight of the fragment is not stable, and thus the cross-sectional area of the 
fragment does not remain constant. To account for this instability, a new cross-sectional 
area is calculated during each pass of the loop. This is accomplished by generating a 
random number using a Normal distribution in which the original cross-sectional area acts 
as the mean, and the standard deviation is 0.05 meters 2 . One might speculate that the drag 
coefficient, C, changes along with the changing cross-sectional area. However, that case 
is not address in this thesis. One thing to remember is that this procedure is simply a 
means to simulate the unstable flight, and does not represent the actual way in which a 
fragment travels. Upon completion of a fragment's flight path the loop is exited and the 
model conducts damage assessment by locating the fragment's hit position within a target 
subregion. 

4. Target Designation 

Defining a specific target in the model, e.g., a particular populated city in the 
middle east, is too restrictive for purposes of this thesis. The objective is to estimate the 
damage to some designated target region caused by fragments generated from an 
intercepted TBM. The model must be adaptable for different targets. To accomplish this, 
however, one restriction is necessary. The dimensions of the target region must be 
rectangular. This rectangular region is then divided into equal-sized blocks, or subregions, 
to allow finer detail in estimating target damage. The number of blocks depends on the 
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desired number of rows and columns requested from the user at the start of the simulation. 
The generic target region used in this example is depicted in Figure 2. The exact 
center of the rectangle is designated as the aim point. This rectangle then acts as a 
damage template which is placed over the target region with the bottom of the rectangle, 
blocks 13 through 16, being perpendicular to the TBM flight path. Figure 3 displays a 
target region composed of blocks 6, 7, 10 and 11 If the intended target is not rectangular 
in shape or is comprised of a number of smaller targets, i.e., runways at a coastal airport 
or piers at a naval port, the size of the rectangle is increased to encompass the total target 
area and only the blocks containing the actual targets are considered areas of interest. 
Figure 4 depicts piers at a naval port composed of blocks 5, 6, 7, 9, 11, 13, and 15. The 
remaining blocks are considered areas of insignificance. Flowever, a target that is located 
within a template block is required to occupy the entire block, otherwise the damage 
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Figure 2. Layout of generic target region template 
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assessment is inaccurate because the damage suffered is not uniform throughout a block. 
This encompassing feature of the target rectangle allows the target to be of any size. 
However, smaller blocks produce more accurate damage estimations because of the finer 
detail. Assuming that all targets are rectangle is fairly unrealistic, but using this 
assumption for damage assessment only, allows this kind of model to be used for just 
about every target imaginable that fits within the outlined target region. 

5. Damage Assessment 

This model assesses damage by three measures: the mean number of fragments to 
impact the target subregions, the mean total mass to impact the target subregions, and the 
mean total kinetic energy to impact the target subregions. The final positions of each 
fragment generated from a single intercept are compared against the location of the 
designated target region through a series of if-then statements to determine the extent of 
damage A realistic way of determining damage with this damage template is to assume 
that the entire block is a target, i.e., there are no insignificant portions of the target block. 
Since this model deals with a generic target region, a generalization is made about the 
damage that occurs to a target block. If a target occupies an entire block then the damage 
that occurs to that block is treated as one amount even though fragments may impact at 
several different points within the block. If it is determined that the fragment impacted the 
target region, the estimated damage totals of the target block which contains that 
fragment's final position are updated. The counter for the number of fragments to impact 
is incremented, the accumulator for the total mass damage has the mass of the fragment 
added to it, and the accumulator for the total kinetic energy damage has the kinetic energy 
of the fragment added to it. The kinetic energy associated with fragment i is calculated by 
the following: 

KE, = \niiv] (20) 
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where 


v = velocity of fragment / at impact 
m, = mass of fragment / at impact 


In order to produce mean damage values with small standard errors, at least 2000 
replications of the model are necessary At the completion of each replication the totals 
for each block are recorded to be used in calculating the means and the standard errors of 
the means using the following formulas: 


R 


D lb = IS £>»(/) 

K J =1 


( 21 ) 


2 = 7T-T £(£>»« - D ®) 2 

K- 1 J= \ 


( 22 ) 



where 


(23) 


Djb(j ) = type i damage value for block b in the jth simulated replication 
D it, = mean type i damage value for block b 
R = number of replications 
s = sample variance of damage value 
SE = standard error of the mean 


In addition to estimating the extent of damage to the target region, the model also 
calculates the amount of damage that occurs to the target's surrounding areas. This thesis 
considers these areas, "miss" areas. However, analysis is conducted in order to know the 
location of the remaining fragment damage in case less significant targets are located in 
these areas. Figure 5 displays these areas. Since the launch point could be from any 
direction, assigning the surrounding areas with a direction, i.e., north of the target, is not 
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implemented. Instead, the model determines if the damage is past or over the target, or is 
short of it, or right or left of it 
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Past & Right 

Left 

TARGET 

Right 

Short & Left 

Short 

Short & Right 


Figure 5. Plot of surrounding areas to the target region. 


6. Structure of a Test 

Many assumptions were made in this model concerning the fragmentation of a 
TBM upon interception. The reason for this is the lack of available data on both the true 
fragmentation and the characteristics of the resulting fragments. One possible way of 
validating this model is to structure a test that replicates an intercepted TBM and the 
resulting fragments. After completion of such a test, a distribution fitted to the fragment 
impact points determines if the assumptions made were correct leading to a valid model. 
An example of such a test, on a small and inexpensive scale, could use 2.75" rockets with 
ball bearings and a radio-controlled burster warhead as test vehicles. At predetermined 
points along the rocket's trajectory the radio-controlled burster is detonated releasing the 
ball bearings to travel their own paths towards the target region, or in this case the witness 
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area. This area is covered in saran wrap enabling exact determination of the ball bearing 
impact points. The positions of these impact points are then fitted to a probability 
distribution to see if the model assumptions are correct. This is just one example of a 
possible test that can be used to validate the proposed model in this thesis since exact data 
is not available. 

B. DESIGN OF A SIMULATION 

A series of runs are conducted to estimate the damage that occurs to subregions 
within a designated target region from fragments generated by the intercept of a TBM as 
the two launch angles, vertical and off-target, are varied slightly to account for the 
different impact points associated with the TBM. To increase the coverage of this model 
and examine other aspects of the damage estimates, three different intercept ranges are 
simulated. Varying the intercept range allows an understanding of the need to intercept 
the TBM as soon as possible, thereby causing less damage to the target area. The 
MODSIM source code for the simulation is contained in Appendix C. 
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IV. ILLUSTRATIVE EXAMPLE 


To aid in the comprehension of this model an illustrative example is provided. The 
TBM is launched towards a target located 560 km downrange, The initial conditions for 
the TBM described earlier are used: vertical launch angle of 30.0 degrees, off-target 
launch angle of 0.0 degrees, launch acceleration of 0.023875 meters/seconds 2 , cross- 
sectional area of 0.6082 meters 2 , boost phase time interval of 90 seconds, and missile 
mass of 6000 kg after boost phase. These parameters achieve a direct hit of the aim point 
if the TBM is not intercepted. The model output is highly sensitive to changes in these 
numbers. 

The additional information necessary to conduct this simulation are the intercept 
range measured in kilometers downrange from the aim point, an expected number of 
fragments generated, and the length and width of the target measured in kilometers. For 
this example, the intercept ranges of interest from the aim point are 2, 3, and 4 kilometers, 
the expected number of fragments is 5 [Ref. 8], and 2 km for the length and width of the 
target, corresponding to a square target area of 4 km 2 . The target region is divided into 
16 equal blocks, or square subregions, each of area 0.25 km 2 . Table 1 lists the estimates 
of the expected number of fragments to impact each block, the expected total mass to 
impact each block and the expected kinetic energy to impact each block, along with the 
standard errors associated with each of the estimates for an intercept of 2 km. Plots of 
these results are displayed in Figures 6 through 11. Table 2 show the results for an 
intercept of 3 km, and the plots of these results are shown in Figures 12 through 17. 

Table 3 shows the results for an intercept of 4 km, and the plots of these results are 
displayed in Figures 18 through 23. The remaining vertical launch angles examined in this 
thesis are the set (29.8°, 29.9°, 30.1°, 30.2°, 30.3°). The remaining off-target launch 
angles are the set ( 0.025°, 0.05°, 0.075°, 0.10°, 0.15°). The off-target launch angles 
represent either side of the target line. The model was run for the right side of the target 
line. To get the data for angles to the left just transpose the data to the corresponding 
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blocks. All possible combinations of these angles were analyzed at the three intercepts 
and the results are contained in Appendix B 



Table 1. List of estimates of expected values and corresponding standard errors of 
the estimates for number of fragment impacts, total mass impact and total 
kinetic energy impact to the target area resulting from a 2 km intercept. 
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Figure 6. Mean number of fragments to impact the target which is divided into 16 
equal blocks and its surrounding areas. This data is generated using 
launch angles of 30.0 degrees vertical and 0.0 degrees off-target and 
an intercept range of 2 km. 


X= Downrange Distance 
Y= Off-Target Distance 



Figure 7. 3-D plot of mean number of fragment impacts per target block utilizing 
launch angles of 30.0 degrees vertical and 0.0 degrees off-target and 
an intercept range of 2 km. 
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- .are 8. Mean total mass (kg) to impact the target divided into 16 equal blocks and 
its surrounding areas. This data is generated using launch angles of 30.0 
degrees vertical and 0.0 degrees off-target and an intercept range of 2 km. 


X= Downrange Distance 
Y= Off-Target Distance 



1258.15 kg 


Figure 9. 3-D plot of mean total mass (kg) impact per target block utilizing launch 
angles of 30.0 degrees vertical and 0.0 degrees off-target and an intercept 
range of 2 km. 
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Figure 10. Mean total kinetic energy (J) to impact the target divided into 16 equal 
blocks and its surrounding areas. This data is generated using launch 
angles of 30.0 degrees vertical and 0.0 degrees off-target and an 
intercept range of 2 km. 
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Figure 11. 3-D plot of mean kinetic energy (J) impact per target block utilizing 

launch angles of 30.0 degrees vertical and 0.0 degrees off-target and an 
intercept range of 2 km. 
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Block Number Expected Value Standard Error 

10-Frag 

0.273 

0.007 

10-Mass 

822.03 kg 

22.894 kg 

10-KE 

12.182 J 

0.394 J 

11-Frag 

0.273 

0.007 

11-Mass 

822.03 kg 

22.894 kg 

ll-KE 

12.182 J 

0.394 J 

■m 

0.471 

0.011 

14-Mass 

712.543 kg 

16.144 kg 

14-KE 

5.749 J 

0.134 J 

15-Frag 

0.471 

0.011 

15-Mass 

712.543 kg 

16.144 kg 

15-KE 

5.749 J 

0.134 J 

Short-Frag 

3.137 

0.061 



29.438 kg 

Short-K E 

510.411 J 

11.848 J 


Table 2. List of estimates of expected values and corresponding standard errors of 
the estimates for number of fragment impacts, total mass impact and total 
kinetic energy impact to the target area resulting from a 3 km intercept. 
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Figure 12. Mean number of fragments to impact the target which is divided into 16 
equal blocks and its surrounding areas. This data is generated using 
launch angles of 30.0 degrees vertical and 0.0 degrees off-target and 
an intercept range of 3 km. 


X= Downrange Distance 



Figure 13. 3-D plot of mean number of fragment impacts per target block utilizing 
launch angles of 30.0 degrees vertical and 0.0 degrees off-target and 
an intercept range of 3 km. 
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Figure 14. Mean total mass (kg) to impact the target divided into 16 equal blocks and 
its surrounding areas. This data is generated using launch angles of 30.0 
degrees vertical and 0.0 degrees off-target and an intercept range of 3 km. 


X= Downrange Distance 


Y= Off-Target Distance 


Z= Mean Mass Impact 



Figure 15. 3-D plot of mean total mass (kg) impact per target block utilizing launch 
angles of 30.0 degrees vertical and 0.0 degrees off-target and an intercept 
range of 3 km. 
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Figure 16. Mean total kinetic energy (J) to impact the target divided into 16 equal 
blocks and its surrounding areas. This data is generated using launch 
angles of 30.0 degrees vertical and 0.0 degrees off-target and an 
intercept range of 3 km. 
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Figure 17. 3-D plot of mean kinetic energy (J) impact per target block utilizing 

launch angles of 30.0 degrees vertical and 0.0 degrees off-target and an 
intercept range of 3 km. 
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Block Number Expected Value Standard Error 

10-Frag 

0.017 

0.002 

1 O’Mass 

88.549 kg 

10.683 kg 

10-KE 

1.435 J 

0.174 J 

11-Frag 

0.017 

0.002 

11 -Mass 

88.549 kg 

10.683 kg 

11-KE 

1.435 J 

0.174 J 

14-Frag 

0.259 

0.007 

14-Mass 

740.46 kg 

21.234 kg 

14-KE 

7.956 J 

0.241 J 

15-Frag 

0.259 

0.007 

15-Mass 

740.46 kg 

21.234 kg 

15-KE 

7.956 J 

0.241 J 

Short-Frag 

4.073 

0.065 

Short-Mass 

2928.78 kg 

42.802 kg 

Short-K E 

460.854 J 

10.895 J 


Table 3. List of estimates of expected values and corresponding standard errors of 
the estimates for number of fragment impacts, total mass impact and total 
kinetic energy impact to the target area resulting from a 4 km intercept. 
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Figure 18. Mean number of fragments to impact the target which is divided into 
16 equal blocks and its surrounding areas. This data is generated using 
launch angles of 30.0 degrees vertical and 0.0 degrees off-target and 
an intercept range of 4 km. 


X— Downrange Distance 


Y= Off-Target Distance 
Z— Mean Number of Hits 



4.073 


Figure 19. 3-D plot of mean number of fragment impacts per target block utilizing 
launch angles of 30.0 degrees vertical and 0.0 degrees off-target and 
an intercept range of 4 km. 
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Figure 20. Mean total mass (kg) to impact the target divided into 16 equal blocks 


and its surrounding areas. This data is generated using launch angles of 
30.0 degrees vertical and 0.0 degrees off-target and an intercept range 
of 4 km. 
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Figure 21. 3-D plot of mean total mass (kg) impact per target block utilizing 

launch angles of 30.0 degrees vertical and 0.0 degrees off-target and 
an intercept range of 4 km. 
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Figure 22. Mean total kinetic energy (J) to impact the target divided into 16 equal 
blocks and its surrounding areas. This data is generated using launch 
angles of 30.0 degrees vertical and 0.0 degrees off-target and an 
intercept range of 4 km. 
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Figure 23. 3-D plot of mean kinetic energy (J) impact per target block utilizing launch 
angles of 30.0 degrees vertical and 0.0 degrees off-target and an 
intercept range of 4 km. 
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To aid in data interpretation, it is assumed that the target damage template 
encompasses the entire target region, meaning each block is a block of interest. Also, 
damage calculations were conducted for the surrounding areas in order to distinguish 
where the remaining fragments fall. An interesting result is that the total number of 
fragment impacts, including the surrounding areas, experienced in the three intercept 
ranges, as shown in Figures 4, 6, and 10, are equal. This is expected since 2000 
replications are run and all impacts are tabulated. The total masses are also the same for 
the three intercept ranges. However, as expected, the total kinetic energies decrease as 
intercept range increases because the terminal velocities of the fragments are smaller. 

Also of note is the symmetry experienced around the target region's centerline when the 
off-target launch angle of 0.0° is used. This result is expected when taking a large enough 
sample. However, with only 2000 replications a small degree of sampling error remained. 
To account for this error, the damage values that resulted for the blocks on either side of 
the target line were averaged together and this result was applied to both blocks. This 
procedure is justified because of the symmetry; it would not be necessary if larger samples 
were taken but without doing so the computing time required would be increased 
unnecessarily. Any other off-target launch angle does not produce symmetry around the 
centerline as shown in Figure 24. The total number of fragment impacts is still the same 
with this new off-target launch angle, but the fragment impact positions have shifted over 
to the right, as anticipated. Increasing this angle more causes a further shift to the right of 
the fragment impact positions up to a point where no fragments impact the target region. 

Determining the actual resulting damage if the TBM was left alone to impact the 
target region at the aim point is very complicated. Knowledge of any structures contained 
within the region would be necessary, along with the manner in which the blast from the 
warhead affects these structures. To simplify this problem, the same damage estimates 
used for the fragments are employed. Because four blocks surround the aim point when 
launch angles of 30.0° vertical and 0.0° off-target are used, the resulting damage is spread 
equally among the four blocks as shown in Figures 25 and 26. 
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Figure 24. Mean number of fragments to impact the target which is divided into 16 
equal blocks and its surrounding areas. This data is generated using 
launch angles of 30.0 degrees vertical and 0.025 degrees ofT-target and 
an intercept range of 2 km. 
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Figure 25. Mean total mass to impact target region when TBM is not 
intercepted and aim point is impacted. 
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Figure 26. Mean total kinetic energy to impact target region when TBM is not 
intercepted and aim point is impacted 
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V. DISCUSSION AND CONCLUSIONS 


The estimates of damage occurring to the target region vary depending on three 
factors The first is the vertical launch angle. As this angle is increased the number of 
fragments that impact the target region increases, along with the respective masses, since 
the TBM is at a higher altitude at intercept than when a smaller angle is employed. This 
is shown in Figures 6 and 27. The higher altitude produces a narrower dispersion pattern 
which concentrates the fragments impacts closer to the target. However, the kinetic 
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Figure 27. Mean number of fragments to impact the target which is divided into 16 
equal blocks and its surrounding areas. This data is generated using 
launch angles of 30.1 degrees vertical and 0.0 degrees off-target and 
an intercept range of 2 km. 

energies are smaller because the fragments fly longer than when intercepted at a lower 
altitude thus allowing air resistance more time to reduce the velocities. The second factor 
is the off-target launch angle. As this angle is increased, i.e., moved to the right of target 
line, the extent of damage to the target region shifts to the right up to a point where all the 
fragments miss the target region. A shift to the left occurs when this angle is decreased, 
i.e., moved to the left of the target line. Since the off-target component of velocity is 
treated equally in both cases, moving right or left of the target line, the extent of damage is 
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equal, however, the block which suffers the damage is transposed to the opposite side of 
the target line. As shown in Appendix B, 0.15° is the angle which shifts everything outside 
of the target region. The shift is shown in Figures 6 and 24. The last factor which causes 
variance in the damage estimates is the intercept range. Intercepts which occur close to 
the aim point generate more fragment impacts to the target region than intercepts which 
are not as close. The difference is illustrated in Figures 6 and 18. 

These three factors, taken independently, yield fairly distinct distributions of 
damage. However, when examined simultaneously, these factors create a very diverse 
damage distribution. Appendix B contains the layouts of all possible combinations of 
these three factors for the example given in the last chapter. The only exceptions are the 
off-target angles to the left of the target line. The extent of damage is the same, but the 
blocks in which the damage occurs are transposed on the opposite side of the target line. 
One possible approach to estimating the extent of damage when these three factors are 
known is to specify a response function, 

Expected Damage = f(0 0 T,®v,R ) (26) 


where 

&OT = off-target launch angle 
0 v = vertical launch angle 

R = intercept range measured in km downrange from the target 

Such a function, specified on physical grounds, could be fitted to data obtained at various 
angles and ranges, and then used to estimate responses at others. 

The conclusion is that knowing the expected damage to occur to a target region if 
a TBM is intercepted is helpful in making a decision whether to intercept it or not. This is 
obvious, but the model is developed to illustrate the necessity of this information. This 
model is a simple analytical attempt to produce a tentative estimate. At this time, very 
little is known concerning the dispersion of fragments generated by an intercept. Because 
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of this lack of information, a crude estimate is the only estimate available. When the 
assumptions made in this thesis are determined with certainty, a more detailed version of 
the model, replicating the exact flight paths and impact points, can be developed. 

Another area that may be of interest is determining which type of intercept 
approach is better, a blast-fragment intercept or a kinetic energy intercept. A blast- 
fragment intercept such as the one simulated in this thesis does not require the accuracy of 
a kinetic energy intercept where the interceptor must actually impact the TBM. However, 
the kinetic energy intercept almost completely destroys the TBM in flight resulting in very 
little chance of collateral damage. The technology required for a kinetic energy 
interceptor is more costly than that of the blast-fragment interceptor, and because of this 
high cost it does not presently exist. The Navy is experimenting with the Light 
Exo-Atmospheric Projectile (LEAP), a kinetic kill vehicle that impacts a TBM outside of 
the atmosphere. If this experiment is successful the technology may be applied to work 
within the atmosphere on intercepts such as the ones proposed in this thesis eliminating 
the need for an estimate of fragment damage. 
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APPENDIX A. 


CALCULATION OF TBM FLIGHT PATH 

As mentioned in the text, only four of the actual characteristics of the A1 Hussein 
are used in this model: the TBM's mass, the TBM's length, the TBM's cross-sectional 
area, and the TBM's boost-phase time interval. The mass of the TBM at launch is 7000 
kg. However, at the end of boost phase the mass is reduced to 6000 kg due to expended 
fuel. The method for compensating for expended fuel is not entirely realistic since the 
mass does not just disappear all at once. However, since the boost phase is not of concern 
in this thesis, and since the assumption has been made that straight-line motion during 
boost was unaffected by air resistance and gravitation, the manner in which missile mass is 
decreased during fuel use is insignificant so long as the mass at intercept equals the mass 
at the completion of boost phase. The length of the TBM is 12 meters. The diameter of 
the TBM is 0.88 meters at its widest point yielding a cross-sectional area of 0.6082 
meters 2 . Finally, determining the exact boost-phase time interval is extremely difficult and 
requires numerous factors, especially the fuel burn rate of the TBM which depends heavily 
upon the density and temperature of air. Usually, as is the case with the A1 Hussein, the 
length of boost phase time interval is given by an approximation. Because of this 
complexity in determining the exact time interval, the approximate boost phase time 
interval of 90 seconds is used in this model [Ref. 2], 

To determine the remaining parameters, an assumption about the distance to the 
target region is required; it is that the center of the target region designated for attack, i.e., 
the aim point, is located 560 km from the launch point. Having knowledge of the target 
distance and choosing an arbitrary launch angle of 30 degrees enables the launch 
acceleration to be calculated with the assumption that acceleration was constant 
throughout the entire boost phase. The calculated launch acceleration is 0.023875 
meters/seconds 2 . The known parameters, together with these parameters calculated 
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under the above assumptions, are the only data required to model the launch of the TBM 
and begin the simulation. Since these calculated characteristics are not the precise 
characteristics of the A1 Hussein TBM, the flight paths modeled in this simulation do not 
replicate those of an AJ Hussein exactly, and thus, are considered generic TBM flight 
paths usefully close to the actual to provide insight. All of the above numbers can be 
changed in the software at the will of an analyst. 

TBM Flight Path Calculations 

Prior to commencing the simulation, the user is prompted for three items that aid 
in the construction of the simulated TBM flight path. The first item is the range at which 
intercept is going to occur, measured in kilometers downrange from the aim point. The 
last two items, the vertical and the off-target launch angles, determine the location of the 
TBM's impact point if it is not intercepted. In this thesis only one pair of angles, 30 
degrees vertical and 0 degrees off-target, result in a direct hit of the aim point. All other 
pairs yield different impact points that are off-target-region center. Launching a TBM, 
even one with state-of-the-art technology, at a target 560 km away results in some level of 
uncertainty as to the exact impact point. Of course, better technology leads to smaller 
uncertainty. Varying the vertical and off-target launch angles slightly, i.e., up to +/- 0.3 
degrees, allows parametric sensitivity analysis to be conducted in this model. Such 
numbers are illustrative and can be altered by the analyst. 

The TBM flight is broken up into two segments. The first segment, which consists 
of the boost phase, is fairly easy to simulate taking into account the assumption of 
straight-line motion, during which air resistance and gravity have negligible effect. Basic 
equations of physics are employed to determine the TBM's position and velocity at the 
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completion of boost phase. The linear distance traveled by the TBM during boost phase, 
D, is determined from a side calculation where gravity and thrust are taken into account. 
From this the three-dimensional position of the TBM at the completion of the boost phase 
is calculated using the following equations: 

X = D cos (LA ) sin (OA ) (1) 

Y = D sin (LA) (2) 

Z = D cos (LA ) cos (OA ) (3) 


where 


X = downrange displacement from launch point 
Y = altitude 

Z = off-target displacement 
LA = vertical launch angle 
OA = off-target angle 

The TBM's velocity in the direction tangent to flight path after boost phase, V 0 , is also 
determined from a side calculation where gravity and thrust are taken into account. From 
this the components of velocity are calculated in the following manner 

Vx - Vq cos (LA) cos (OA) (4) 


Vy = Vo sin ( LA ) 
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Vz = Vo cos (LA ) sin (OA ) (6) 

where 

V x = downrange component of velocity 
V v = vertical component of velocity 
V z = off-target component of velocity 

The second segment consists of the TBM ballistic flight path. The purpose for its 
generation is to enable examination of various intercept points so as to determine the 
lethality of the fragments created by those intercepts to parts of the target region. 
Therefore, it is necessary to be able to calculate the positions and velocities of the TBM 
throughout its entire flight path. Initially, an attempt is made to replicate the exact ballistic 
path the TBM would follow through the atmosphere. Knowing that gravity and air 
resistance are the only two forces which affect a projectile's flight, an expression of air 
resistance is sought to implement in a closed-form differential equation of motion. It has 
been determined that retardation of a projectile due to air resistance is a function of the 
projectile's velocity [Ref. 3], 

The actual flight of the TBM is modeled using three-dimensional equations of 
motion taking gravity and atmospheric density into consideration. In the simulation, a 
loop is created which continues until the downrange displacement of the TBM equals the 
intercept range entered by the user. The position and velocity of the TBM, calculated by 
Equations (1) through (6), are the initial conditions prior to entering the loop. For ease of 
computation, an interval of one second is chosen as the time step in the loop. At each 


48 



APPENDIX A. 


CALCULATION OF TBM FLIGHT PATH 


time step the TBM's terminal velocity at time t is calculated by 


V T {t) = 


Cp(tM 


where 


m = mass of the TBM 

g = acceleration due to gravity 

C = drag coefficient dependent on angle of attack 

p(t) = density of air which is a function of altitude at time t 

A = cross-sectional area of the TBM 


Since the remaining equations are all functions of time, the reference to time, i.e., (t), is a 
given and thus removed from the equations. The velocity of the TBM at each pass 
through the loop is updated using 


V = (V 2 X + V\ + v 2 z )i 


From these two calculations, the changes in each component of velocity are calculated 
using the following equations of motion 
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£V 

dt 


g 


VV\ 


Y g 


( 10 ) 


cil 


W 7 .& 


( 11 ) 


From these equations much simpler forms for finding the velocity components at each 
(1-second) time step are defined by 


Vx« 


Vy„ 




Wx t CpA 

~ ^- 


nycpA 


__ VV^CyA 
2m 


( 12 ) 

(13) 

(14) 


The simple one-step differential equation solver portrayed above can be replaced by far 
more accurate numerical procedures, such as 4th-order Runge-Kutta, or other 
methodologies in standard package programs [Ref. 4], In view of the exploratory nature 
of this investigation such complexity has not been introduced. 

Atmospheric density, p, is extremely difficult to model. The Earth's atmosphere is 
primarily composed of nitrogen and oxygen. Solar radiation affects the dynamic 
properties of this medium by constantly changing the temperature, pressure, and chemical 
constituents, particulate presence and electrical properties [Ref. 5], The inability to model 
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atmospheric density creates a problem when trying to replicate the exact flight path of the 
TBM. Therefore, the values for p utilized in this model are taken from the U.S. Standard 
Atmosphere Table [Ref. 6], The drag coefficient, C, takes into account the relative 
contributions of viscous and form resistances, and it depends on the nature of the TBM 
(size, shape, and irregularity and roughness of its surface) as well as on the characteristics 
of the flow of air over the TBM [Ref. 7], Numerous factors are involved in computing C; 
a value of 1.0 is used in the present illustrative simulation. 

The position of the TBM is also required at every step in the loop. Having 
calculated the velocity, the position can be determined rather easily by 


Xi 

ii 

> 

+ 

(15) 

Y t 

II 

i 

> 

+ 

(16) 

z, 

ii 

N 

7 

> 

+ 

(17) 


The loop continues until the downrange displacement, X,, equals the intercept point 
determined by the user. When the intercept point is reached the TBM separates into a 
random number of fragments. 
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RESULTS OF SIMULATION WITH VARIOUS LAUNCH ANGLES 
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RESULTS OF SIMULATION WITH VARIOUS LAUNCH ANGLES 
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Intercept Range = 4.0 km 


Each Block Contains: # of Fragments 
Total Mass 
Total KE 










































APPENDIX B. 

RESULTS OF SIMULATION WITH VARIOUS LAUNCH ANGLES 



Vertical Launch Angle = 30.0° Each Block Contains: # of Fragments 


Off-Target Launch Angle = 0.10° 


Total Mass 


Intercept Range = 4.0 km 


Total KE 



Vertical Launch Angle = 30.0° Each Block Contains: # of Fragments 

Off-Target Launch Angle = 0.15° Total Mass 


Intercept Range = 4.0 km 


Total KE 


97 


























APPENDIX B. 


RESULTS OF SIMULATION WITH VARIOUS LAUNCH ANGLES 


P&L 

0 


Past 

0 



P& R 

0 


0 

0 

0 

0 


Left 

0 

0.020 

100.606 

1.385 

0.020 

100.606 

1.385 

0 

Right 

0 

0 

0.177 

551.304 

5.686 

0.177 

551.304 

5.686 

0 

0 


0 

0.312 

582.163 

4.678 

0.312 

582.163 

4.678 

0 


S&L 

0 


Short 

2.609 

2118.654 

634.174 



S&R 

0 


Vertical Launch Angle = 30.1° Each Block Contains: # of Fragments 

Off-Target Launch Angle = 0.0° Total Mass 

Intercept Range = 4.0 km Total KE 


P&L 

0 


Past 

0 



P&R 

0 


0 

0 

0 

0 


Left 

0 

0 

0.040 

201.212 

2.770 

0 

Right 

0 

0 

0 

0.355 

1102.608 

11.372 

0 

0 


0 

0 

0.623 

1164.325 

9.355 

0 


S&L 

0 


Short 

3.609 

2118.654 

634.174 



S&R 

0 


Vertical Launch Angle = 30.1° 
Off-Target Launch Angle = 0.025° 
Intercept Range = 4.0 km 


Each Block Contains: # of Fragments 
Total Mass 
Total KE 


98 














APPENDIX B. 


RESULTS OF SIMULATION WITH VARIOUS LAUNCH ANGLES 


P&L 

; o 


Past 

0 



P&R 

0 


0 

0 

0 

0 


Left 

0 

0 

0.028 

140.454 

1.936 

0.012 

60.757 

0.834 

Right 

0 

0 

0 

0.241 

745.640 

7.670 

0.114 

356.967 

3.703 

0 


0 

0 

0.456 

843.933 

6.756 

0.168 

320.392 

2.600 


S & L 

0 


Short 

3.609 

2118.654 

634.174 



S&R 

0 

Vertical Launch Angle = 30.1° 
Off-Target Launch Angle = 0.05° 
Intercept Range = 4.0 km 

Each Block Contains: # of F 
Total 
Total 

P&L 

0 


Past 

0 



P&R 

0 


0 

0 

0 

0 


Left 

0 

0 

0 

0.040 

201.212 

2.770 

Right 

0 

0 

0 

0 

0.354 

1101.383 

11.362 

0 


0 

0 

0 

0.624 

1165.550 

9.366 


S&L 

0 


Short 

3.609 

2118.654 

634.174 



S&R 

0 


Vertical Launch Angle = 30.1° 
Off-Target Launch Angle = 0.075° 
Intercept Range = 4.0 km 


Each Block Contains: # of Fragments 
Total Mass 
Total KE 


99 




















APPENDIX B. 


RESULTS OF SIMULATION WITH VARIOUS LAUNCH ANGLES 



Vertical Launch Angle = 30.1° Each Block Contains: # of Fragments 

Off-Target Launch Angle = 0.10° Total Mass 


Intercept Range = 4.0 km 


Total KE 



Vertical Launch Angle = 30.1° Each Block Contains: # of Fragments 

Off-Target Launch Angle = 0.15° Total Mass 

Intercept Range = 4.0 km Total KE 


100 

























APPENDIX B. 


RESULTS OF SIMULATION WITH VARIOUS LAUNCH ANGLES 


P & L 


Past 



P&R 

0 


0 



0 


0 

0 

0 

0 




0.012 

0.012 




0 

66.104 

66.104 

0 


Left 


0.812 

0.812 


Right 

0 


0.098 

0.098 


0 


0 

355.413 

355.413 

0 




3.655 

3.655 





0.228 

0.228 




0 

529.662 

529.662 

0 




4.567 

4.567 



S&L 


Short 



S&R 

0 


3.950 



0 



2684.442 






471.931 





Vertical Launch Angle = 30.2° Each Block Contains: # of Fragments 

Off-Target Launch Angle = 0.0° Total Mass 


Intercept Range = 4.0 km 


Total KE 


P&L 

0 


Past 

0 



P&R 

0 


0 

0 

0 

0 


Left 

0 

0 

0.025 

132.207 

1.624 

0 

Right 

0 

0 

0 

0.196 

710.826 

7.310 

0 

0 


0 

0 

0.455 

1058.373 

9.126 

0 


S&L 

0 


Short 

3.951 

2685.392 

471.939 



S&R 

0 


Vertical Launch Angle = 30.2° 
Off-Target Launch Angle = 0.025° 
Intercept Range = 4.0 km 


Each Block Contains: # of Fragments 
Total Mass 
Total KE 


101 




















APPENDIX B. 


RESULTS OF SIMULATION WITH VARIOUS LAUNCH ANGLES 
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APPENDIX C. 


MODSIM SIMULATION SOURCE CODE 


DEFINITION MODULE NewFrag; 

FROM NewTBM IMPORT G. PositionType; 

CONST C= 1.0; 

TYPE TargetType - RECORD 
Rows: INTEGER: 

Columns: INTEGER: 

END RECORD: 

IntArrayType = ARRAY INTEGER. INTEGER OF INTEGER; 

RealArrayType = ARRAY INTEGER. INTEGER OF REAL; 

FragObj = OBJECT 
Mass : REAL; 

Area : REAL; 

OriginalArea: REAL; 

Velocity : PositionType; 

Position : PositionType; 

CurrentVel: REAL; 

Damage : REAL; 

LengthofBlock: REAL; 

Southwest: PositionType; 

DX: REAL; 

DZ: REAL; 

TGT: TargetType; 

ASK METHOD Objlnit; 

ASK METHOD GetMass(IN M: REAL); 

ASK METHOD CalcArea; 

ASK METHOD CalcVelocity(IN V: PositionType); 

ASK METHOD GetPosition(IN P: PositionType); 

ASK METHOD SetTargetRegion(IN LE: REAL; IN WI: REAL; IN RO: INTEGER; IN CO: INTEGER; IN DI: 

REAL); 

ASK METHOD CalcDamage; 

ASK METHOD Travel; 

END OBJECT; 

FragmentsType = ARRAY INTEGER OF FragObj; 

VAR 

Fragment: FragObj; 

Fragments : FragmentsType; 

FragCount: IntArrayType; 

MassTotal : RealArrayType; 

KETotal : RealArrayType; 

END MODULE. 
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MODSIM SIMULATION SOURCE CODE 


DEFINITION MODULE NewRunRcp; 
PROCEDURE RunReplications: 

END MODULE. 


DEFINITION MODULE NewTBM; 

FROM RandMod IMPORT RandomObj; 

FROM NcwFrag IMPORT IntArrayType, RealArrayType; 

CONST G = .0098; 

Accci = .023875; 

Area - 6082; 

WM = 7000.0; 

MBT - 90.0; 
seedl = 1; 
seed2 = 3; 
seed3 = 5; 
seed4 = 7; 

Cl = 1.0; 

TYPE PositionTvpe = RECORD 
X; REAL; 

Y: REAL; 

Z; REAL; 

END RECORD; 

MissileObj = OBJECT 
Dist; REAL; 

WarheadMass; REAL; 

LA; REAL; 

OA: REAL; 

TBMPosition: PositionType; 

TBMVelocity: PositionType; 

Acceleration; REAL; 

Pieces; INTEGER; 

MotorBurnTime: REAL; 

MissileArea; REAL; 

InterceptRange; REAL; 

AfterMotorSpeed; REAL; 

InterceptTime: REAL; 

Altitude: REAL; 

MeanFragments: REAL; 

Length: REAL; 

Width: REAL; 

Rows: INTEGER; 
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MODSIM SIMULATION SOURCE CODE 


Columns: INTEGER: 

TotalFragCount: IntArrayTvpe; 

TotalMD: ReaiArrayTvpe; 

TotalKE; ReaiArrayTvpe: 

TPLD: REAL; 

TPLMD: REAL; 

TPLKE: REAL; 

TSLD: REAL; 

TSLMD: REAL; 

TSLKE: REAL; 

TLD: REAL; 

TLMD: REAL; 

TLKE: REAL; 

TPRD: REAL; 

TPRMD: REAL; 

TPRKE: REAL; 

TSRD: REAL; 

TSRMD: REAL; 

TSRKE: REAL; 

TRD: REAL; 

TRMD: REAL; 

TRKE: REAL; 

TPD: REAL; 

TPMD; REAL; 

TPKE: REAL; 

TSD: REAL; 

TSMD: REAL; 

TSKE: REAL; 

ASK METHOD Objlnit; 

ASK METHOD FlyTraj; 

ASK METHOD Breakup: 

ASK METHOD GetInterceptRange(IN R; REAL); 

ASK METHOD GetMeanFragment(IN ME: REAL); 

ASK METHOD GetAngles(IN VA; REAL; IN OTA: REAL); 

ASK METHOD GetAltitude(IN ALT: REAL); 

ASK METHOD GetTargetDimension(IN L: REAL; IN W: REAL; IN R: INTEGER; IN C: INTEGER; IN D: 
REAL); ASK METHOD ResetDamage(IN Zero: REAL); 

ASK METHOD GetDamage(IN FC: INTEGER; IN MT: REAL; IN KET: REAL; IN i: INTEGER; IN j: 
INTEGER); ASK METHOD GetPastLeftDamage(IN DD:REAL; IN DMD:REAL; IN DKE:REAL); 

ASK METHOD GetShortLeftDamage(IN DD:REAL; IN DMD:REAL; IN DKE:REAL); 

ASK METHOD GetLeftDamage(IN DD:REAL; IN DMD:REAL; IN DKE:REAL); 

ASK METHOD GetPastRightDamage(IN DD:REAL; IN DMD:REAL; IN DKE:REAL); 

ASK METHOD GetShortRightDamage(IN DD:REAL; IN DMD:REAL; IN DKE:REAL); 

ASK METHOD GetRightDamage(IN DD:REAL; IN DMD:REAL; IN DKE:REAL); 

ASK METHOD GetPastDamage(IN DD:REAL; IN DMD:REAL; IN DKE:REAL); 

ASK METHOD GetShortDamage(IN DD:REAL; IN DMD:REAL; IN DKE:REAL); 
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MODSIM SIMULATION SOURCE CODE 


END OBJECT. 

VAR TBM: MissilcObj: 
P.R.S.T: RandomObj; 
END MODULE. 


DEFINITION MODULE NewWriteLine: 
PROCEDURE WriteLine(IN Siring: STRING): 
END MODULE. 


IMPLEMENTATION MODULE NewFrag; 

FROM RandMod IMPORT RandomObj; 

FROM NewTBM IMPORT PositionType, G, TBM.R,S,T; 
FROM MathMod IMPORT SQRT.POWER. pi,CEIL; 
FROM NewWriteLine IMPORT WriteLine; 


OBJECT FragObj: 


ASK METHOD Objlnit: 

i 

\ . 

BEGIN 

NEW(Posilion); 
NEW(Velocity); 
NEW(TGT);* 
NEW(SouthWest); 

END METHOD: 


} 

} 


! . } 

ASK METHOD GetMass(IN M: REAL); 

..} 

BEGIN 

Mass := M; 

END METHOD: 


a 

i. 

ASK METHOD CalcArea; 

{ .. 

VAR 

Radius: REAL: 


} 

} 
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MODSIM SIMULATION SOURCE CODE 


BEGIN 

Radius := ((3.0/(4.0*pi)) * (0.001216425 * Mass)); 
Area := pi * POWER(Radius,(2.0/3.0)); 
OriginaiArea := Area; 

END METHOD; 


ASK METHOD CalcVeIocitv(IN V: PositionTvpe); 

{ . } 

VAR 

kickx.kickv.kickz ; REAL; 

BEGIN 

kickx ;= (ASK T TO Normal(O.O.O.Ol)); 
kicky := (ASK T TO Normal(0.0.0.01)); 
kickz := (ASK T TO Normal(0.0 r 0.01)); 

Velocity.X := V.X + kickx; 

Velocity.Y := V.Y + kicky; 

Velocity.Z := V.Z + kickz; 

END METHOD; 


ASK METHOD GetPosition(IN P: PositionType); 

i .} 

BEGIN 

Position.X ;= P.X; 

Position. Y := P.Y; 

Position.Z := P.Z; 

END METHOD; 


ASK METHOD SetTargetRegion(IN LE; REAL; IN WI: REAL; IN RO; INTEGER; IN CO; INTEGER; IN DI: 

REAL); 

{ . > 

BEGIN 

DX := LE/FLOAT(RO); 

DZ := WI/FLOAT(CO); 

SouthWest.X := DI - (0.5*FLOAT(RO)*DX); 

SouthWest.Z := 0.0 - (0.5*FLOAT(CO)*DZ); 

TGT.Rows := RO; 

TGT.Columns := CO: 

NEW(FragCount. 1..RO. 1..CO); 

NEW(MassTotal. 1..RO, 1..CO); 

NEW(KETotal. 1..RO, 1..CO); 

END METHOD; 
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MODSIM SIMULATION SOURCE CODE 


i I 

l. f 

ASK METHOD CalcDamage: 

{ . } 

VAR 

IJ. PLD.PD.PRD,SLD.SD,SRD.LD.RD: INTEGER; 

PLMD.PLKE.PMD.PKE.PRMD.PRKE ? SLMD,SLKE.SMD.SKE. SRMD,SRKE ? LMDXKE,RMD,RKE: REAL; 
BEGIN 

I := CElL((Position.X - South'West.X)/DX); 

J '= CEIL((Position.Z - SouthWcst.Z)/DZ); 

IF (I >0) AND (I <= TGT.Rows) 

IF (J > 0) AND (J <= TGT.Columns) 

INC(FragCount[IJ]); 

MassTotalJIJ] := MassTotal[IJ] + Mass; 

KETotal[I ? J] ~ KETotal[IJ] + (0.5*Mass*POWER(CurrentVel f 2.0)); 

ASK TBM TO GetDamage(FragCount[I,J],MassTotal[I,J],KETotal[I,J],I,J); 

END IF; 

END IF; 

IF I > TGT.Rows 
IF J <= 0 
PLD := 1; 

PLMD := Mass; 

PLKE ;= (0.5*Mass*POWER(CurrentVeL2.0)); 

ASK TBM TO GetPastLeftDamage(FLOAT(PLD),PLMD ? PLKE); 

END IF; 

IF (J > 0) AND (J<= TGT.Columns) 

PD := 1; 

PMD := Mass; 

PKE := (0.5*Mass*POWER(CurrentVeL2.0)); 

ASK TBM TO GetPastDamage(FLOAT(PD),PMDEKE); 

END IF; 

IF J > TGT.Columns 
PRD := 1; 

PRMD := Mass; 

PRKE ;= (0.5*Mass*POWER(CurrentVel f 2.0)); 

ASK TBM TO GetPastRightDamage(FLOAT(PRD),PRMD,PRKE); 

END IF; 

END IF; 

IF I <= 0 
IF J <= 0 
SLD := 1; 

SLMD := Mass; 

SLKE := (0.5*Mass*POWER(CurrentVel,2.0)); 

ASK TBM TO GetShortLeftDamage(FLOAT(SLD),SLMD,SLKE); 

END IF; 

IF (J > 0) AND (J <= TGT.Columns) 

SD := 1; 

SMD := Mass; 
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SKE := (0.5*iMass*POWER(CurrentVel.2.0)); 

ASK TBM TO GetShortDamage(FLOATfSD) ? SMD,SKE); 

END IF; 

IF J>TGT.Columns 

SRD := 1; 

SRMD := Mass; 

SRKE := (0.5*Mass*POWER(CurrentVeF2.0)); 

ASK TBM TO GetShortRightDamage(FLOAT(SRD),SRMD,SRKE); 
END IF; 

END IF; 

IF (I > 0) AND (I <= TGT.Rows) 

IF J <= 0 
LD := 1; 

LMD := Mass; 

LKE ;= (0.5*Mass*POWER(CurrentVel,2.0)); 

ASK TBM TO GetLeftDamage(FLOAT(LD).LMD,LKE); 

END IF; 

IF J >= TGT.Columns 
RD := 1; 

RMD := Mass; 

RKE := (0.5*Mass*POWER(CurrentVel,2.0)); 

ASK TBM TO GetRightDamage(FLOAT(RD),RMD,RKE); 

END IF; 

END IF; 

END METHOD; 

( 

i. 

ASK METHOD Travel 


VAR 

t. rho.VTerm : REAL; 

BEGIN 
t := 0.0; 

WHILE (Position.Y > 0.0) AND (Velocity.X > 0.0) 
CASE ROUND(Position.Y) 

WHEN 0 .1: 
rho := 1.1117; 

WHEN 2; 
rho ;= 1.0066; 

WHEN 3: 
rho := 0.90925; 

WHEN 4: 
rho := 0.81935; 

WHEN 5: 
rho := 0.73643; 

WHEN 6; 
rho := 0.66011; 
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WHEN 7: 
rho := 0.59002; 
WHEN 8: 
rho := 0.52579; 
WHEN 9: 
rho := 0.46706; 
WHEN 10: 
rho := 0.41351; 
WHEN II; 
rho := 0.3648; 

WHEN 12; 
rho := 0.31194; 
WHEN 13: 
rho := 0.2666; 

WHEN 14: 
rho := 0.22786; 
WHEN 15: 
rho := 0.19475; 
WHEN 16: 
rho := 0.16647; 
WHEN 17: 
rho := 0.1423; 

WHEN 18: 
rho := 0.12165; 
WHEN 19: 
rho := 0.104; 

WHEN 20: 
rho := 0.08991; 
WHEN 21..25: 
rho := 0.04; 

WHEN 26..30: 
rho := 0.01841; 
WHEN 31..40: 
rho := 0.0039957; 
WHEN 41..50: 
rho := 0.0010269; 
WHEN 51..60: 
rho := 0.00030592; 
WHEN 61..70: 
rho := 0.000087535; 
WHEN 71..80: 
rho :- 0.00001999; 
WHEN 81..150: 
rho := 0.0000004974; 
END CASE; 
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CurrentVel := SQRT(POWER(Velocity.X.2.0) + POWER(Velocity.Y,2.0) + POWER(Velocity.Z,2.0)); 
Velocity.X := ((Velocity.X* 1000.0) - ((Current Vel* 1000.0*Velocity.X*l()00.0*rho*Area*C)/(2.0*Mass))) 
/ 1000 . 0 ; 

Velocity. Y := ((Velocity. Y*1000.0) - (G*1000.0) - ((Current Vel* 1000.0* Velocity. Y*1000.0*rho*Area*C)/ 
(2.0*Mass)))/1000.0; 

Velocity.Z := ((Velocity.Z* 1000.0) - ((CurrentVel* 1000.0*Velocity.Z*1000.0*rho*Area*C)/(2.0*Mass))) 
/ 1000 . 0 : 

Position.X := Position.X + (Velocity.X); 

Position.Y := Position.Y + (Velocity.Y); 

Position.Z := Position.Z + (Velocity.Z); 

Area := ASK S TO Normal(OriginalArea. 0.05); 
t :=t+ 1.0; 

END WHILE; 

ASK SELF TO CalcDamage; 

END METHOD; 

END OBJECT; 

END MODULE. 


IMPLEMENTATION MODULE NewRunRep; 

FROM NewTBM IMPORT TBM; 

FROM SimMod IMPORT StartSimulation,ResetSimTime; 

FROM NewWriteLine IMPORT WriteLine; 

FROM MathMod IMPORT POWER, SQRT; 

FROM NewFrag IMPORT RealArrayType; 

PROCEDURE RunReplications; 

VAR 

Z.p.q.c.d.ROW.COL: INTEGER; 

APLD.APLMD.APLKE.ASLD. ASLMD,ASLKE,ALD.ALMD,ALKE,APRD,APRMD,APRKE,ASRD: REAL; 
ASRMD.ASRKE.ARD.ARMD,ARKE,APD,APMD,APKE,ASD,ASMD,ASKE,PLD2,PLMD2,PLKE2: REAL; 
SLD2.SLMD2 SLKE2,LD2,LMD2,LKE2,PRD2,PRMD2,PRKE2,SRD2,SRMD2,SRKE2,RD2; REAL; 
RMD2.RKE2.PD2.PMD2.PKE2,SD2,SMD2,SKE2: REAL; 
K,T,Range,VertAngle,OTAngle,MeanFragment,Len,Wid,Distance: REAL; 
AvgFragCount.AvgMass,AvgKE,SQFragCount,SQMass,SQKE: RealArrayType; 

BEGIN 

NEW(TBM); 

UTPUT("What is the downrange intercept range from the target?"); 

INPUT(Range); 

OUTPUT("What is the desired vertical launch angle?"); 

INPUT(VertAngle); 

OUTPUT("What is the desired off-target launch angle?"); 

INPUT(OTAngle); 

OUTPUT("What is the mean number of fragments desired?"); 

INPUT(MeanFragment); 
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OUTPUT("What is the length of the target in KM?"); 
lNPUT(Len); 

OUTPUT("What is the width of the target in KM?"); 

INPUT(Wid); 

OUTPUT("How many rows of blocks in the target region?"); 

INPUT(ROW); 

OUTPUT("How many columns of blocks in the target region?"); 

INPUT(COL); 

OUTPUTC What is the distance to the aimpoint in the center of the target?"); 

INPUT(Distance); 

ASK TBM TO GetMeanFragment(MeanFragntent); 

ASK TBM TO GetAngles(VertAngle.OTAngle); 

ASK TBM TO GclTargetDimension(Len.Wid.ROW.COL.Distance); 

ASK TBM TO GetlnterccptRange(Range); 

NEW(AvgFragCount. 1..ROW. 1..COL); 

NEW(AvgMass. 1..ROW. 1..COL); 

NEW(AvgKE. 1. ROW. 1..COL); 

N E W( SQF ragCount. 1 ..ROW. 1 ..COL); 

NEW(SQMass. 1..ROW. 1..COL); 

NEW(SQKE. 1..ROW. 1..COL); 

Z := 1; 

WHILE Z< 101 
ASK TBM TO FlvTraj; 

StartSimulation; 

FOR c ;= 1 TO ROW 
FOR d := 1 TO COL 

AvgFragCount[c.d] := AvgFragCount[c,d] + FLOAT(ASK TBM TotalFragCount[c,d]); 
AvgMass[c.d] ;= AvgMass[c,d] + (ASK TBM TotalMD[c,d]); 

AvgKEMl := AvgKEfc.d] + (ASK TBM TotalKE[c,d]); 

SQFragCountlc.dJ := SQFragCount[c,dJ + POWER(FLOAT(ASK TBM TotalFragCount[c,d]),2.0); 
SQMass[c,d] := SQMass[c.d] + POWER(ASK TBM TotalMD[c.d],2.0); 

SQKEfc.d] := SQKE[c.d] + POWER(ASK TBM TotalKE[c,d],2.0); 

END FOR; 

END FOR; 

APLD := APLD + (ASK TBM TPLD); 

APLMD := APLMD + (ASK TBM TPLMD); 

APLKE := APLKE + (ASK TBM TPLKE); 

ASLD := ASLD + (ASK TBM TSLD); 

ASLMD := ASLMD + (ASK TBM TSLMD); 

ASLKE := ASLKE + (ASK TBM TSLKE); 

ALD := ALD + (ASK TBM TLD); 

ALMD := ALMD + (ASK TBM TLMD); 

ALKE := ALKE + (ASK TBM TLKE); 

APRD := APRD + (ASK TBM TPRD); 

APRMD := APRMD + (ASK TBM TPRMD); 

AJPRKE := APRKE + (ASK TBM TPRKE); 

ASRD := ASRD + (ASK TBM TSRD); 

ASRMD := ASRMD + (ASK TBM TSRMD); 
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ASRKE := ASRKE + (ASK TBM TSRKE); 

ARD := ARD + (ASK TBM TRD); 

ARMD := ARMD + (ASK TBM TRMD); 

ARKE := ARKE + (ASK TBM TRKE); 

APD := APD + (ASK TBM TPD): 

APMD := APMD + (ASK TBM TPMD): 

APKE := .APKE + (ASK TBM TPKE); 

ASD := ASD + (ASK TBM TSD); 

ASMD := ASMD + (ASK TBM TSMD); 

ASKE := ASKE + (ASK TBM TSKE); 

PLD2 := PLD2 + POWER(ASK TBM TPLD.2.0); 

PLMD2 := PLMD2 + POWER(ASK TBM TPLMD.2.0); 

PLKE2 := PLKE2 + POWER(ASK TBM TPLKE.2.0); 

SLD2 := SLD2 + POWER(ASK TBM TSLD.2.0); 

SLMD2 := SLMD2 + POWER(ASK TBM TSLMD.2.0); 

SLKE2 := SLKE2 + POWER(ASK TBM TSLKE,2.0); 

LD2 := LD2 + POWER(ASK TBM TLD,2.0); 

LMD2 := LMD2 + POWER(ASK TBM TLMD.2.0); 

LKE2 := LKE2 + POWER(ASK TBM TLKE,2.0); 

PRD2 = PRD2 + POWER(ASK TBM TPRD.2.0); 

PRMD2 := PRMD2 + POWER(ASK TBM TPRMD.2.0); 

PRKE2 := PRKE2 + POWER(ASK TBM TPRKE.2.0); 

SRD2 := SRD2 + POWER(ASK TBM TSRD,2.0); 

SRMD2 := SRMD2 + POWER(ASK TBM TSRMD,2.0); 

SRKE2 := SRKE2 + POWER(ASK TBM TSRKE,2.0); 

RD2 := RD2 + POWER(ASK TBM TRD,2.0); 

RMD2 := RMD2 + POWER(ASK TBM TRMD,2.0); 

RKE2 := RKE2 + POWER(ASK TBM TRKE,2.0); 

PD2 := PD2 + POWER(ASK TBM TPD.2.0); 

PMD2 := PMD2 + POWER(ASK TBM TPMD.2.0); 

PKE2 := PKE2 + POWER(ASK TBM TPKE,2.0); 

SD2 := SD2 + POWER(ASK TBM TSD,2.0); 

SMD2 := SMD2 + POWER(ASK TBM TSMD,2.0); 

SKE2 := SKE2 + POWER(ASK TBM TSKE,2.0); 

ASK TBM TO ResetDamage(O.O); 

ResetSimTime(O.O); Z := Z + 1; 

END WHILE; 

K := FLOAT(Z-l); 

T := K - 1.0: 

WriteLinef "); 

WriteLine("Intercept Distance is: " + REALTOSTR(Range) + " km"); 
WriteLine("Vertical Launch Angle is: " + REALTOSTR(VertAngle) + " degrees"); 
WriteLine("Off-Target Launch Angle is: " + REALTOSTR(OTAngle) + " degrees"); 
WriteLineC'Intercept Altitude is: " + REALTOSTR(ASK TBM Altitude) + " km"); 
WriteLinef "): 
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FOR p = i TO ROW 
FOR q = 1 TO COL 

AvgFragCount(p.q] := AvgFragCount[p.q]/K; 

AvgMass|p.q| := AvgMass(p.q]/K: 

AvgKJE(p.q) := AvgKE|p.qJ/K: 

SQFragCountip.q] := SQRT(((SQFragCount[p,q]/T) - ((K/T)*POWER(AvgFragCount[p,q],2.0)))/K); 

SQMass|p.q] := SQRT(((SQMass[p.q]/T) - ((K/T)*POWER(AvgMass[p,q],2.0)))/K); 

SQKE|p.q] := SQRT(((SQKE[p.q]/T) - ((K/T)*POWER(AvgKE[p,q],2.0)))/K); 

WritcLinc("Block " + INTTOSTR(p) + + INTTOSTR(q) + " MEAN HITS: " + 

REALTOSTR(AvgFragCount|p.q]) + " Std Error: " + REALTOSTR(SQFragCount[p.q])); 

WriteLinc("Block " + INTTOSTR(p) + + INTTOSTR(q) + " MEAN MASS: " + REALTOSTR(AvgMass[p,q]) 

- "Std Error: " + REALTOSTR(SQMass(p.q])); 

WriicLmc("Block " - INTTOSTR(p) + + INTTOSTR(q) + " MEAN KE: " + REALTOSTR(AvgKE[p,qj) 

- " Std Error: " + REALTOSTR(SQKE[p.q])); 

END FOR: 

END FOR: 


APLD := APLD/K; 

APLMD := APLMD/K: 

APLKE := APLKE/K; 

ASLD := ASLD/K; 

ASLMD := ASLMD/K: 

ASLKE := ASLKE/K; 

ALD := ALD/K; 

ALMD := ALMD/K: 

ALKE := ALKE/K; 

APRD := APRD/K: 

APRMD := APRMD/K: 

APRKE := APRKE/K: 

ASRD := ASRD/K; 

ASRMD := ASRMD/K: 

ASRKE := ASRKE/K; 

ARD := ARD/K; 

ARMD := ARMD/K: 

ARJKE := ARKE/K: 

APD := APD/K; 

APMD := APMD/K: 

APKE := APKE/K: 

ASD := ASD/K; 

ASMD := ASMD/K; 

ASKE := ASKE/K; 

PLD2 := SQRT((PLD2/T)-((K/T)*POWER(APLD,2.0)))/SQRT(K); 
PLMD2 := SQRT((PLMD2/T)-((K/T)*POWER(APLMD,2.0)))/SQRT(K); 
PLKE2 := SQRT((PLKE2/T)-((K/T)*POWER(APLKE,2.0)))/SQRT(K); 
SLD2 := SQRT((SLD2/T)-((K/T)*POWER(ASLD,2.0)))/SQRT(K); 
SLMD2 := SQRT((SLMD2/T)-((K/T)*POWER(ASLMD,2.0)))/SQRT(K); 
SLKE2 := SQRT((SLKE2/T)-((K/T)*POWER(ASLKE,2.0)))/SQRT(K); 
LD2 := SQRT((LD2/T)-((K/T)*POWER(ALD,2.0)))/SQRT(K); 
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LMD2 := SQRT((LMD2/T)-((K/T)*POWER(ALMD.2.0)))/SQRT(K); 

LKE2 := SQRT((LKE2/T)-((K/T)*POWER(ALKE.2.0)))/SQRT(K); 

PRD2 := SQRT((PRD2/T)-((K/T)*POWER(APRD.2.0)))/SQRT(K); 

PRMD2 := SQRT((PRMD2/T)-((K/T)*POWER(APRMD.2.0)))/SQRT(K); 

PRKE2 := SQRT((PRBCE2/T)-((BC/T)*POWER(APRBCE,2.0)))/SQRT(K); 

SRD2 := SQRT((SRD2/T)-((K/T)*POWER(ASRD,2.0)))/SQRT(K); 

SRMD2 := SQRT((SRMD2/T)-((K7T)*POWER(ASRMD,2.0)))/SQRT(K); 

SRKE2 := SQRT((SRKE2/T)-((K/T)*POWER(ASRKE.2.0)))/SQRT(K); 

RD2 := SQRT((RD2/T)-((BC/T)*POWER(ARD.2.0)))/SQRT(K); 

RMD2 := SQRT((RMD2/T)-((K/T)*POWER(ARMD,2.0)))/SQRT(K); 

RKE2 := SQRT((RKE2/T)-((K/T)*POWER(ARKE.2.0)))/SQRT(K); 

PD2 := SQRT((PD2/T)-((BC/T)*POWER(APD,2.0)))/SQRT(K); 

PMD2 := SQRT((PMD2/T)-((BC/T)*POWER(APMD.2.0)))/SQRT(K); 

PBCE2 := SQRT((PKE2/T)-((K/T)*POWER(APKE.2.0)))/SQRT(K); 

SD2 := SQRT((SD2/T)-((K/T)*POWER(ASD.2.0)))/SQRT(K); 

SMD2 := SQRT((SMD2/T)-((K/T)*POWER(ASMD.2.0)))/SQRT(K); 

SKE2 := SQRT((SBCE2/T)-((K/T)*POWER(ASBCE.2.0)))/SQRT(K); 

WriteLine("MISSES/LOCATIONS"); 

WriteLineC"'); 

WriteLine("Mean PastLeft Misses: " + REALTOSTR(APLD) + " Std Error: " + REALTOSTR(PLD2)); 
WriteLine("Mean PastLeft Mass Damage: " + REALTOSTR(APLMD) + " Std Error: " + 

REALTOSTR(PLMD2)); 

WriteLinefMean PastLeft KE Damage: " + REALTOSTR(APLKE) + " Std Error: " + REALTOSTR(PLKE2)); 
WriteLinefMean ShortLeft Misses: " + REALTOSTR(ASLD) + " Std Error: " + REALTOSTR(SLD2)); 
WriteLineC'Mean ShortLeft Mass Damage:" + REALTOSTR(ASLMD) + " Std Error: "+REALTOSTR(SLMD2)); 
WriteLineC'Mean ShortLeft KE Damage: " + REALTOSTR(ASLKE) + " Std Error: "+ REALTOSTR(SLKE2)); 
WriteLineC’Mean Left Misses: " + REALTOSTR(ALD) + " Std Error: " + REALTOSTR(LD2)); 

WriteLineC'Mean Left Mass Damage: " + REALTOSTR(ALMD) + " Std Error: "+ REALTOSTR(LMD2)); 
WriteLineC'Mean Left KE Damage: " + REALTOSTR(ALKE) + " Std Error: " + REALTOSTR(LKE2»; 

WriteLineC'Mean PastRight Misses: " + REALTOSTR(APRD) + " Std Error: " + REALTOSTR(PRD2)); 
WriteLineC'Mean PastRight Mass Damage:" + REALTOSTR(APRMD) + " Std Error: "+ 

REALTOSTR(PRMD2)); WriteLineC'Mean PastRight KE Damage: " + REALTOSTR(APRKE) + " Std Error: " 

+ REALTOSTR(PRKE2)); WriteLineC'Mean ShortRight Misses: " + REALTOSTR(ASRD) + " Std Error: " + 
REALTOSTR(SRD2)); WriteLineC'Mean ShortRight Mass Damage:" + REALTOSTR(ASRMD) + " Std Error: 
"+REALTOSTR(SRMD2)); WriteLineC'Mean SkortRight KE Damage: " + REALTOSTR(ASRKE) + " Std Error: 
" + REALTOSTR(SRKE2)); WriteLineC'Mean Right Misses: " + REALTOSTR(ARD) + " Std Error:" + 

REALTOSTR(RD2)); WriteLineC'Mean Right Mass Damage: " + REALTOSTR(ARMD) + " Std Error: " + 
REALTOSTR(RMD2)); WriteLineC'Mean Right BCE Damage: " + REALTOSTR(ARBCE) + " Std Error:" + 
REALTOSTR(RKE2)); WriteLineC'Mean Past Misses: " + REALTOSTR(APD) + " Std Error: " + 

RE ALTOSTR(PD2)); WriteLineC'Mean Past Mass Damage: " + REALTOSTR(APMD) + " Std Error: " + 
REALTOSTR(PMD2)); WriteLineC'Mean Past BCE Damage: " + REALTOSTR(APBCE) + " Std Error: " + 

REALTOSTR(PBCE2)); WriteLineC’Mean Short Misses: " + REALTOSTR(ASD) + " Std Error:" + 

REALTOSTR(SD2)); WriteLineC'Mean Short Mass Damage: " + REALTOSTR(ASMD) + " Std Error: " + 
REALTOSTR(SMD2)); WriteLineC'Mean Short BCE Damage: " + REALTOSTR(ASBCE) + " Std Error: " + 
REALTOSTR(SBCE2)); 
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END PROCEDURE: 
END MODULE. 


IMPLEMENTATION MODULE NewTBM; 

FROM NewWriteLine IMPORT WriteLine: 

FROM SimMod IMPORT SimTime: 

FROM RandMod IMPORT RandomObj.FctchSeed: 

FROM NcwFrag IMPORT FragObj. FragmentsTvpe, IntArravType, RealArravType: 
FROM MathMod IMPORT COS, S IN, SQRT, POWER, pi; 

OBJECT MissileObj; 

{ . } 

ASK METHOD Objlnit: 

{ . } 

BEGIN 

NEW(TBM Velocity); 

NEW(TBMPosition): 

NEW(P): 

NEW(R): 

NEW(S): 

NEW(T): 

ASK P TO SetSeed(FetchSeed(seedl)); 

ASK R TO SetSeed(FetchSeed(seed2)); 

ASK S TO SetSeed(FetchSeed(seed3)); 

ASK T TO SetSeed(FetchSeed(seed4)); 

Acceleration := Accel; {kni/sec2} 

WarheadMass := WM; {kg} 

MissileArea := Area; {m A 2} 

MotorBurnTime := MBT : {secs} 

AfterMotorSpeed := Acceleration * MotorBurnTime; 

END METHOD: 

: . } 

ASK METHOD FlyTraj; 

{ . } 

VAR 

InterceptPt, rho. CurrentVel: REAL; 

BEGIN 

NEW(TotalFragCount. L.Rows, 1..Columns); 

NEW(TotalMD. L.Rows. 1 ..Columns): 

NEW(TotalKE. L.Rows. 1..Columns); 

InterceptPt := Dist - InterceptRange; 
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TBMPosition.X := (Accel*0.5*MBT*MBT*COS(LA)*COS(OA)); 

TBMPosition.Y •= (Accel*0.5*MBT*MBT*SIN(LA)); 

TBMPosition.Z := (Acce!*0.5*MBT*MBT*COS(LA)*SIN(OA)); 

TBMVcIocity.X := (Accel*MBT)*COS(LA)*COS(OA); 

TBMVelocity.Y := (Accel*MBT)*SIN(LA); 

TBMVelocity.Z := (Accel*MBT)*COS(LA)*SIN(OA); 

WarheadMass := 6000.0; 

WHILE FLOAT(ROUND(TBMPosition.X)) < InterceptPt 

CurrentVel := SQRT(POWER(TBMVelocity.X,2.0)+POWER(TBMVelocitv.Y,2.0) + 
POWER(TBM Velocity.Z.2.0)); 

IF TBMPosition.Y > 0.0 
CASE ROUND(TBMPosition.Y) 

WHEN ().. 1; 
rho =1 1117; 

WHEN 2: 
rho .= 1.0066; 

WHEN 3: 
rho := 0.90925; 

WHEN 4; 
rho := 0.81935; 

WHEN 5: 
rho := 0.73643; 

WHEN 6: 
rho := 0.66011; 

WHEN 7: 
rho := 0.59002; 

WHEN 8: 
rho := 0.52579; 

WHEN 9; 
rho := 0.46706; 

WHEN 10: 
rho := 0.41351; 

WHEN 11: 
rho := 0.3648; 

WHEN 12: 
rho := 0.31194; 

WHEN 13: 
rho := 0.2666; 

WHEN 14: 
rho := 0.22786: 

WHEN 15: 
rho := 0.19475; 

WHEN 16: 
rho := 0.16647; 

WHEN 17: 
rho := 0.1423; 
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WHEN 18: 
rho := 0.12165: 

WHEN 19: 
rho := 0.104: 

WHEN 20: 
rho := 0.08991: 

WHEN 21..25: 
rho := 0.04: 

WHEN 26..30: 
rho := 0.01841: 

WHEN 31..40: 
rho := 0.0039957: 

WHEN 41..50: 
rho := 0.0010269: 

WHEN 51..60: 
rho := 0.00030592: 

WHEN 61..70: 
rho := 0.000087535: 

WHEN 71..80: 
rho := 0.00001999: 

WHEN 81..200: 
rho := 0.0000004974: 

END CASE: 

END IF; 

TBMVelocity.X := ((TBMVelocity.X* 1000.0) - ((Current Vel* 1000.0* TBMVelocitv.X* 1000.0* 
rho*MissileArea*Cl) /(2.0*WarheadMass)))/1000.0; 

TBMVelocity. Y := ((TBMVelocity.Y* 1000.0) - (0.0098*1000.0) - ((CurrentVel* 1000.0*TBMVelocity.Y* 
1000.0*rho*MissileArea*Cl)/(2.0*WarheadMass)))/1000.0; 

TBMVelocirv.Z := ((TBMVelocity.Z* 1000.0) - ((CurrentVel* 1000.0*TBMVelocitv.Z* 1000.0* rho* 
MissileArea* Cl)/(2.0*WarheadMass)))/1000.0; 

TBMPosition.X := TBMPosition.X + TBMVelocity.X; 

TBMFosition.Y := TBMPosition.Y + TBM Velocity. Y; 

TBMPosilion.Z := TBMPosition.Z + TBMVelocitv.Z; 

END WHILE: 

ASK SELF TO Breakup; 

ASK SELF TO GetAltitude(TBMPosition. Y); 

END METHOD; 

{ . } 

ASK METHOD Breakup; 

( . } 

VAR 

l. K: INTEGER: 

m. TotalMass: REAL: 

Fragment: FragObj: 

Fragments: FragmentsType; 
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APPENDIX C. 


MODSIM SIMULATION SOURCE CODE 


BEGIN 

TolalMass := 0.0; 

K ■= ASK R TO Poisson(ASK SELF MeanFragments); 

IF K <= I 
K :=2; 

END IF; 

NEW(Fragments, 1..5*K); 

I := 1; 

WHILE TotalMass < WarheadMass 

m := (WarheadMass/FLOAT(K))*(ASK P TO Exponential(l.O)); 

TotalMass := TotalMass -+- m; 

IF m > 70.0 

IF TotalMass < WarheadMass 
NEW(FragmeiuslIJ); 

ASK Fragments) I] TO GetMass(m); 

ASK Fragments!I) TO CalcArea; 

ASK Fragments[I] TO CalcVelocitv(TBMVelocity); 

ASK Fragments[I] TO GetPosition(TBMPosition); 

ASK Fragmentsfl] TO SetTargetRegion(Length,Width,Rows,Columns, Dist); 
ASK Fragments[I] TO Travel; 

I := I + 1; 

END IF; 

END IF; 

END WHILE; 

END METHOD; 


! . } 

ASK METHOD GetlnterceptRange(IN R: REAL); 

{.} 

BEGIN 

nterceptRange := R; 

END METHOD; 


{.} 

ASK METHOD GetMeanFragment(IN MF: REAL); 


> 

* 


BEGIN 


} 


MeanFragments := MF; 
END METHOD; 


{.} 

ASK METHOD GetAngles(IN VA: REAL; IN OTA: REAL); 

{ . } 

BEGIN 

LA := (VA*pi)/180.0; 

OA := (OTA*pi)/180.0; 

END METHOD; 
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APPENDIX C 


MODSIM SIMULATION SOURCE CODE 


< ) 

l./ 

ASK METHOD GeiAltitude(IN ALT: REAL); 

{.} 

BEGIN 

Altitude = ALT: 

END METHOD: 

!.} 

ASK METHOD GclTargetDimension(IN L: REAL: IN W: REAL: IN R: INTEGER; IN C: INTEGER: IN D 

REAL); 

* > 

l.I 

BEGIN 
Length := L; 

Width := W: 

Rows := R; 

Columns := C; 

Dist := D: 

END METHOD: 

;.} 

ASK METHOD GeiDamage(IN FC: INTEGER; IN MT: REAL; IN KET: REAL; IN i: INTEGER; IN j: 

INTEGER); 

..} 

BEGIN 

TotalFragCount[i.j] := TotalFragCountfij] + FC; 

TotalMDfi.j] := TotalMD[ij] + MT; 

TotalKEfij] := TotalKE[ij] + KET; 

END METHOD; 

;.> 

ASK METHOD GetPastLeftDamage(IN DD: REAL; IN DMD: REAL; IN DKE: REAL); 

!.} 

BEGIN 

TPLD ;= TPLD + DD; 

TPLMD := TPLMD + DMD; 

TPLKE := TPLKE + DKE; 

END METHOD; 

(.} 

ASK METHOD GetShortLeftDamage(IN DD: REAL; IN DMD: REAL; IN DKE' REAL)- 

■s.} 

BEGIN 

TSLD := TSLD + DD; 

TSLMD := TSLMD + DMD; 

TSLKE := TSLKE + DKE; 

END METHOD: 
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APPENDIX C. 


MODSIM SIMULATION SOURCE CODE 


ASK METHOD GctLcftDamage(IN DD 


BEGIN 

TLD := TLD + DD; 

TLMD := TLMD + DMD; 

TLKE := TLKE + DKE; 

END METHOD: 

< >. 

t. f 

ASK METHOD GelPastRigluDamage(IN DD: REAL; IN DMD: REAL; IN DKE: REAL); 

i > 

I. f 

BEGIN 

TPRD = TPRD+ DD; 

TPRMD := TPRMD + DMD; 

TPRKE := TPRKE + DKE: 

END METHOD: 

{.} 

ASK METHOD GetShortRightDamage(IN DD: REAL; IN DMD: REAL; IN DKE: REAL); 

{.} 

BEGIN 

TSRD := TSRD + DD: 

TSRMD := TSRMD + DMD; 

TSRKE := TSRKE + DKE; 

END METHOD; 

{.} 

ASK METHOD GetRightDamage(IN DD: REAL; IN DMD: REAL; IN DKE: REAL); 

{ . } 

BEGIN 

TRD := TRD + DD; 

TRMD := TRMD + DMD; 

TRKE := TRKE + DKE; 

END METHOD: 


: REAL; IN DMD: REAL; IN DKE: REAL); 


} 


ASK METHOD GetPastDamage(IN DD: REAL; IN DMD: REAL; IN DKE: REAL); 

{.} 


BEGIN 


TPD := TPD + DD; 
TPMD := TPMD + DMD; 
TPKE := TPKE + DKE; 
END METHOD; 
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APPENDIX C. 


MODSIM SIMULATION SOURCE CODE 


( . } 

ASK METHOD GetShortDamagc(IN DD: REAL; IN DMD: REAL; IN DKE: REAL) 

<. } 

BEGIN 

TSD := TSD + DD; 

TSMD ;= TSMD -r DMD; 

TSKE ;= TSKE *f DKE; 

END METHOD; 


i./ 

ASK METHOD ResetDamage(IN Zero: REAL); 

: . } 

VAR 

k.n: INTEGER; 

BEGIN 

FOR k := 1 TO Rows 
FOR n := 1 TO Columns 
TotalFragCount[k,n] := 0; 

TotalMD[k,n] :=Zero; 

TotalKE[k.n] := Zero; 

END FOR; 

END FOR; 

TPLD := Zero; 

TPLMD := Zero; 

TPLKE := Zero; 

TSLD := Zero; 

TSLMD := Zero; 

TSLKE := Zero; 

TLD := Zero; 

TLMD := Zero; 

TLKE := Zero; 

TPRD := Zero; 

TPRMD := Zero; 

TPRKE := Zero; 

TSRD := Zero; 

TSRMD := Zero; 

TSRKE := Zero; 

TRD := Zero; 

TRMD := Zero; 

TRKE := Zero; 

TPD := Zero; 

TPMD := Zero; 

TPKE := Zero; 

TSD ;= Zero; 

TSMD := Zero; 







APPENDIX C. 


MODSIM SIMULATION SOURCE CODE 


TSKE := Zero; 
END METHOD; 
END OBJECT; 

END MODULE. 


IMPLEMENTATION MODULE NewWriteLine; 

FROM lOMod IMPORT FileUseType(Output); 
FROM IOMod IMPORT StreamObj; 

FROM UtilMod IMPORT DateTime; 

FROM SiniMod IMPORT SimTime; 

FROM NewTBM IMPORT TBM; 

VAR 

DT ; STRING; 

TraceStream ; StreamObj; 

PROCEDURE WriteLine(IN String; STRING); 
BEGIN 

IF (TraceStream = NILOBJ) 

NEW(TraceStream); 

ASK TraceStream TO Open("sim_.out”. Output); 

DateTime(DT); 

ASK TraceStream TO WriteString(DT); 

ASK TraceStream TO WriteLn; 

ASK TraceStream TO WriteLn; E 
ND IF; 

ASK TraceStream TO WriteReal(SimTime(), 7,3); 
ASK TraceStream TO WriteString(" — " + String); 
ASK TraceStream TO WriteLn; 

END PROCEDURE; 

PROCEDURE WriteLineClose; 

BEGIN 

ASK TraceStream TO Close; 

END PROCEDURE; 

END MODULE. 
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